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Abstract 

We review the status of theory and experiment of very rare and forbidden 
kaon decays. We then review the radiative non-leptonic decays, and the as- 
sociated Dalitz pair modes. We pay particular attention to the study of long 
distance physics in radiative decays within the framework of chiral perturbation 
theory (xPT). We discuss the experiments that will run in the near future and 
the modes that they will be able to study. 
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1 Introduction 

Rare kaon decays offer the possibility of probing liigfi energy scales by doing precise 
low energy measurements. In this case, the precision is obtained by looking at very 
rare decays. The advantage of the kaon system is, of course, the long lifetime of both 
the Kl and the K^. In this paper we review the status of theory and experiment of 
rare, forbidden and radiative kaon decays. Other recent reviews include Refs. ^ 
|,|,§,§,0,|,|. 

We concentrate in three distinct classes of decays. We first discuss a set of modes 
that does not occur in the standard model. Studies of these modes constitute searches 
for new physics, and in some cases these rare kaon decays are the most sensitive probes 
for certain kinds of new interactions. 

We then study those modes that occur in the standard model mostly through 
short distance physics, and that are thus amenable to a conventional perturbative 
treatment. These modes occur at the one-loop level, via penguin and box diagrams as 
discussed in section |[ The calculation of these transitions is by now standard[|l^, |lT|, 
and they are dominated by the top-quark intermediate state. The contribution of 
the charm-quark, as well as perturbative QCD corrections, are also known. The 
theoretical aspects of these calculations have been recently reviewed in the literature 
by Buras and Harlander, Ref. 0]. 

The interest of these modes lies in the possibility of measuring some of the pa- 
rameters in the standard model. The rates are sensitive to the values of the top 
quark mass and of the CKM mixing angles. For our discussion, we will assume a 



unitary 3x3 CKM matrix as parameterized by Wolfenstein |T2[. The advantage of 
this approximate parameterization is that it exhibits the hierarchical structure of the 
CKM matrix. It yields expressions that show what is the mixing angle suppression 
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of a transition. This form of the CKM matrix is: 



I-AV2 A AX^p-ir]) 

v= -A 1-AV2 Ax^ n.i] 

AX^{1 - p-iri) -AX^ 1 

The present knowledge of the mixing angles in the CKM matrix, Eq. is summa- 
rized in Ref. g i: 



V 



A 



A 



pi _|_ ^2 



0.22 

0.9 ±0.1 
0.4 ±0.2 

0.59 ±0.18 



;i.2) 

;i.3) 

;i-4) 



Eq. ^]2| is known from Kij, and hyperon decays while Eq. |1.3| is now usually 
extracted from analyses of i? ^ D(D*^ev decays Ijl^. Eq. [L^ is obtained by fitting 
the observed lepton energy spectrum from B semileptonic decay as a sum of 6 — > c 



and 6 — > M contributions [0. The two different numbers result from the new CLEO 
value for Vuh and the from the old CLEO and ARGUS result. To obtain separate 
determinations of p and r], additional input is necessary. However, the other currently 
available sources of information, such as B — B mixing, e and e'/ e are afflicted both 
by dependence on the top-quark mass and calculational difficulties. Since we expect a 
direct measurement of the top-quark mass from collider experiments before too long, 
the latter problem is perhaps more severe. Here rare kaon decays have the potential 
to contribute substantially. The measurement of the decay rates of processes such 
as Ti^uV and Kl tt^uV can provide constraints on p and 77 with very little 

theoretical uncertainty, once the top quark mass is measured. Other measurements 
such as that of the parity-violating decay asymmetry in —>■ n^fi'^fi" can 
provide complementary information, given only modest theoretical input. The first 
measurements of p and 77 will complete our knowledge of the CKM matrix. Subsequent 
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determinations of any of the parameters will permit us to test the three generation 
structure of the standard model. 

Finally we discuss radiative decays that are expected to be dominated by long 
distance physics. In this case we do not have predictions directly from the standard 
model because we do not know how to handle the non-perturbative aspects of the 
strong interactions. Although we expect that this problem will be solved at some 
point by lattice calculations, at present, the only systematic framework we have to 
study these modes is chiral perturbation theory. Within xPT, we parameterize our 
ignorance of strong interaction dynamics in terms of a few unknown coupling con- 
stants. These constants arc then measured in some processes, and after that, they 
can be used to predict additional decay modes. We review the basics of this approach 
in the following section. 
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2 Chiral perturbation theory 

A conventional calculation of non-leptonic kaon decays in the standard model leads 
to the evaluation of matrix elements of four-quark operators between meson states. 
This is a non-perturbative problem that remains unsolved. In xPT one replaces the 
standard model with an effective low energy field theory written directly in terms 
of meson fields. Effective field theories contain an infinite number of operators and 
are, therefore, not very useful unless one has a way to organize the operators and 
to identify the most important ones. xPT organizes the operators in the low energy 
effective Lagrangian in terms of the number of derivatives (and external fields) that 
occur. This corresponds to an expansion of amplitudes in powers of the momentum 
(or energy) of the external particles. The energy scale for this expansion is set by the 
scale of chiral symmetry breaking, Acsb, empirically about 1 GeV O, ll^, 



The effective Lagrangian is constructed on the basis of chiral symmetry, an ap- 
proximate symmetry of the QCD Lagrangian. In the limit of massless u, d and s 
quarks QCD has a global SU{?>)l x SU{?>)r chiral symmetry that is broken spon- 
taneously to SU{'i)v- The chiral Lagrangian is a compact way to keep track of the 
SU{'i)v relations between amplitudes, as well as of the relations between amplitudes 
with different numbers of pions that follow from the spontaneously broken global 
symmetry (the soft pion theorems). 

Apart from including all possible operators that are chirally invariant (organized 
in terms of number of derivatives), xPT incorporates deviations from chiral symmetry 
due to the small quark masses. This leads to other operators that are organized as 
an expansion in powers of meson masses (and number of derivatives). 

Since the typical momentum of particles in a kaon decay is roughly the kaon mass, 
we can give a very rough estimate for the size of the corrections that one can expect 
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in a xPT calculation: they will be of order mj^/l GeV^ ~ 25% of the highest order 
kept in the calculation. 

The framework of xPT has proved extremely useful for analyzing low energy 
processes involving the pseudoscalar meson octet and photons. At low energies, the 
strong and electromagnetic interactions of these particles can be adequately described 
with a chiral Lagrangian with up to four derivatives. The most general chiral La- 



grangian to this order has been written down by Gasser and Leutwyler ||T^. It 
consists of two terms at leading order, C(p^): 

^(2) ^ L^^f^D^UD'^U^^ + fioyTrl^Mf/ + M^. (2.1) 

M is the diagonal matrix {mu,md,ms), and the meson fields are contained in the 
matrix U = exp{2i(f) / f^^) with: 

TT- _7r7V2 + Vv^ ^° (2-2) 

K- 2r]/VQ 
U transforms under the chiral group as t/ — > RUL'^ . We will restrict ourselves to the 
case where photons are the only external fields. In this case the covariant derivative 
is given by (we will not discuss radiative semileptonic decays [pO|] ): 



1 

72 



V 



D^U = d^U -ieA^[Q,U]. 



(2.3) 



and Q is the diagonal matrix (—2/3, 1/3, 1/3). The two constants that appear at this 
order are the pion decay constant pT| , 



U = (92.4 ± 0.2) MeV 



(2.4) 



and the ratio between meson masses and current quark masses Bq. Ignoring isospin 
breaking, = = m. Bo is given by: 
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At this order, the Gell-Mann-Okubo relation follows 

,2 , _2 _ A^2 



3m; + mi = Amj^ (2.6) 



At next to leading order, 0{p^), there are ten more operators [T^ without epsilon 
tensors. For the processes of interest we will only need two out of these ten terms. 
When photons are the only external fields they read: 

C^s^ = -ieL^F^JI^Q (^D^U^D^U + D^UD'^U^^ + eLioF^"" F^^Tr (uQU^Q^ (2.7) 

and F^u is the usual electromagnetic field strength tensor. At this same order there 
are also terms that contain epsilon tensors. These have the same origin as the triangle 
anomaly responsible for vr" 77, and do not involve any unknown coefficients. They 



are contained in the Wess-Zumino-Witten anomalous action |^2|], and the terms of 
interest to us are: 



.(4) _ a p^f 7] 



le 



A^d^n+dPn-d" fyr^ + (2.8) 



A complete calculation to 0{p'^) consists of tree- level diagrams with vertices from 
Eqs. |2.1| , |2.7| , and |2.8| , and of one-loop diagrams using only Eq. |2.1| . The divergences 
that appear in the loop calculation are absorbed by renormalization of the couplings 
in Eq. |2.7| . The renormalized couplings that we will use are defined by regularizing 
in n = 4 — e dimensions. In terms of Aq = (2/e + ln47r + 1 — 7 — ln/i^)/327r^, they are 



T7| : Lg 10 (/u) = -^9,io='='^o/4- These, and the other eight coupling constants have been 



fixed from experiment [T^. This completes the framework needed to discuss strong 



and electromagnetic processes involving the pseudoscalar meson octet and photons. 

To study kaon decays we need, in addition, an effective Lagrangian for the weak 
interactions. In the standard model, the dominant \ AS\ = 1 operators in the effective 
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weak Hamiltonian transform as (8^,, 1r) or {27l, 1r) under chiral rotations. We can 
write a chiral representation for operators with these transformation properties, and 
once again organize them in terms of the number of derivatives. The lowest order 



Lagrangian constructed in this way contains two derivatives p3 

,(2) _ Gf 



r\^} _ _L|T/ T/* 



(78Tr( \L^L^ ) + <7it" ' ( ^Mi3^2i + ^/.23(4L^i + SL^a! 



V2 

+ 5'27 ' M -^^13-^21 + -^m23(-^11 "~ -^22) 



(2.9) 

where = if^UD^W . We can use this Lagrangian to compute K ^ nn decays, 



and fit the unknown constants from experiment [^, ^ The result is well known, 
the amplitudes with |A/| = 3/2 are much smaller than those with |A/| = 1/2. In 
terms of the couplings of Eq. |2]^ the result is: 



+ 5.1 



1^^1=1: ^ 0.16 



32. (2.10) 



15^27 

From now on, we will use \gs\ ~ 5.1 and drop the terms that transform as {27 l, 1/^). 
We will also use the notation: 

Gs = ^\VudV:jgs ^ 9.1 X 10"^ GeV-^ (2.11) 

The situation at next to leading order is much more complicated: a very large 
number of operators, and therefore of unknown coupling constants, has been identified 
||26|| . However, for the radiative decays that we will consider, only a few of those 



operators play a role under the following assumptions ||27| : 



1. Octet dominance. We will not include |A/| = 3/2 operators. 

2. CPS symmetry. The effective weak Hamiltonian in the standard model is invariant 

under a CP transformation followed by the interchange of d ^ s. This same 
symmetry is imposed on the effective Lagrangian. 
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3. Photons are the only external fields, so that F^^, = = eQF^y. 

4. For the normal (odd) intrinsic parity terms (those without (with) an epsilon ten- 

sor), we will be interested only in terms with at most two (three) meson fields. 



The next to leading order weak Lagrangian then reads [^, p8[] : 



J IT 



+ W2TtIqL^X,lJ +ieftwiF^M\QUQU^ 



(2.12) 



for the normal intrinsic parity sector, and 

aiTrf^QL^jTrJ^AeL, 
+ a2TAUQU^L^Ti(\L,\+a^Ti{\[UQU\L^L,] 



J TT 



(2.13) 

for the odd intrinsic parity sector. 

As in the case of the strong interaction Lagrangian, a calculation to 0{p^) involves 
tree and one-loop graphs with Eq. |2.9| , and tree graphs with Eqs. |2.12| , p.l3| . The loop 
graphs are again divergent, and the divergences are absorbed by renormalization of 
the couplings in Eq. |2.12| ||27|| : 2(/^) = 'U^i,2=fAo, wKfi) = w^ — ^Xq. After using a few 



experimental results as input to fit the unknown coefficients, we can proceed to make 
predictions for other processes. A detailed fit of these constants from experiment is 
in itself interesting, as it provides a compact parameterization of low energy data 
that can be used to compare with first principles calculations when these become 
available. In the meantime, they also provide a framework for systematic tests of 
different models. 

There have been some attempts to derive the effective Lagrangian. The case of 
Eq. |2]^ has been studied in detail using a 1/Nc analysis of the strong interactions 
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pl| , and a quark model |3^. There has also been considerable activity using reso- 
nance saturation models and quark models |^ to estimate the 0{p'^) coupling 
constants. We summarize our knowledge of the couplings in the strong interaction 
sector in Table D p5| . 



Table 1: Values of Lg^ 



10- 





Llifi = mp) 


L^o(/i = mp) 


Experiment 


(6.9 ±0.2) X 10-3 


(-5.2 ±0.3) X 10-3 


Vector Dominance 


7.3 X 10-3 


-5.8 X 10-3 


Quark Models 


6.3 X 10-3 


-3.2 X 10-3 



The vector dominance model is a tree-level resonance saturation model. It pro- 
duces scale independent couplings, and the resonance saturation assumption is im- 
plemented by identifying them with the running couplings at a scale equal to the 
resonance mass. There are many variations of the quark model results, but they all 
start from considering a free quark loop. These models are inspired by the 1/Nc 
expansion, but they are not complete calculations to leading order in l/Nc. They 
also result in couplings that are scale independent. We will identify them with the 
running couplings at a scale equal to twice the constituent quark mass, roughly the 
rho mass (the scale dependence is non- leading in l/N^ so we are free to choose any 
scale). Attempts to incorporate gluonic corrections seem to improve the agreement 
with experiment, although precise quantitative predictions are not available |l34 |. 

For the weak interaction parameters there are also several models. Among them 



the weak deformation model l^g, and quark models []37| , p8| . We summarize these 
results in Table 



Similarly, the large-A^^c factorizable contributions to the constants in Eq. |2.13| are 
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Table 2: Model ca 


culations of Wi 










weak deformation model [SG] 


0.007 


0.028 


-0.021 


factorization |37| 


0.025 


0.025 


0.0 


quark model 


-0.003 


0.013 


-0.005 



33: 



1 



ai 



47r2' 87r2' IGtt^' 
For some of the modes that we will discuss, we will need the matrix elements of 
currents; in analogy with the semileptonic decays. For these we will use the following 
form of the currents in terms of mesons: 



03 



(2.14) 



for the charged current, and 



(2.15) 



—I 



(2.16) 



for the neutral current. These results follow from Eq. ^71] [|T^], and are valid to lowest 
order in xPT. We have only kept terms with up to two mesons. From Eq. |2.16| we 
can also find expressions for the scalar and pseudoscalar densities: 

1 



sd 



v2 



S75 



(2.17) 



Bq is the same as that in Eq. |2.5| . Although Bq is not a physical quantity, we will 
simply use Bq/iti^ ^ 15. We have divided by the muon mass for convenience. This 
number reflects that matrix elements of scalar operators are somewhat enhanced. 
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3 Lepton family number violating decays 

In the minimal standard model with massless neutrinos, the lepton family number is 
absolutely conserved so these decays do not occur. The observation of Kl — > yU^e^ 
oi K ^ TTfie decays would therefore be a clear indication of new physics. 

A model independent study of this type of processes can be done following the 



approach of Buchmiiller and Wyler, Ref. |39|. The physics beyond the standard model 



is parameterized by an effective Lagrangian that is gauge invariant under SU{3)c x 
SU{2)i X U{1)y- The leading effective Lagrangian, relevant for our discussion, is 
given by the sum of four fermion operators of the form: 



^ gl- (1 + 75) _ ^ (1 + 75) ^ 

Cv-A^^^l^. es-i^" d 



c/^ _ (1-75) _ (1 + 75) , .on 
Os±p = Cs±pj^fi es d (3.1) 

We have included a (weak) gauge coupling g'^ to reflect the fact that we think of 
these operators as originating in the exchange of a heavy gauge boson (or perhaps a 
scalar) in the new physics sector. We will take Qn = g2 for simplicity and absorb any 
difference, as well as any mixing angles or other factors into the coefficient Cj. The 
factors of 2 have been included for later convenience. In the absence of any dynamical 
information to the contrary, it is natural to assume that Ci is of order 0{1). A is 
the scale that characterizes the heavy degrees of freedom, typically the mass of the 
exchanged boson. It is then conventional to take Cj = 1 and interpret the bounds on 
the decays induced by these operators as bounds on the "scale of new physics" A. We 
have not listed all possible Dirac structures in Eq. pTTI , but just two illustrative ones. 
The operators could be purely vector or axial-vector operators, for example. Tensor 
operators, however, can be reduced to the others by Fierz rearrangements . 



We compare the rates induced by these operators with Standard Model rates using 
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the S.M. effective four-fermion operator for semileptonic decays: 



r c/i . . _ (1 + 75) _ (1 + 75) .0 o^ 

= 2^ "'''^'^ 2 ^^^^ 2 ^ ^ 



and using Eqs. |2l5| , |2l6| , |l7. 



Models that violate lepton flavor number will also induce processes like /i — »• 67, 
/i^ — >■ e^e^e", /i — e conversion in the field of a heavy nucleus, and AM{Ki — Ks) 
among others. For the latter one finds, AM = /l^mj^/SA^ |3^. The experimental 
value AM{Kl-Ks) = (3.522 ±0.016) x lO'^^ j^eV [0, places a bound A > 830 TeV 
for (y — A) ^ {V — A) operators. This bound will be better than the one obtained 
from Kl — >■ fi^e^ until this process reaches a sensitivity around 10~^^. However, this 
comparison of different processes assumes that all the coefficients Ci in the effective 
Lagrangian are of the same order. Although this is a natural assumption, it may not 
be true for given coefficients in specific models. Therefore, it is important to study all 
the different processes, since at some level they provide complementary information. 

3.1 Kl^ fi^e^ 

It is standard to compare this mode to the rate for —>■ fi^u, since in the limit of 
vanishing electron mass it has the same kinematics. Given the pseudoscalar nature of 
the kaon, only the axial vector quark current or pseudoscalar density can contribute, 
as can be seen from Eqs. |2.16| , p.lTj One finds: 



r(K+^/i+z/) \Vus?\g2 
ioT V — A operators, and 



mm 



for scalar, pseudoscalar operators 
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Fig. |1| a and b show the recent results in the search for ^e. The lack of 

a signal in the region corresponding to M(/ie) ~ with aligned initial and final 
states allows 90% c.l. upper limits to be set. These are B{Kl Ai^e^) < 9.4 x 10"^^ 
for KEK-137 and B{Kl fi^e^) < 3.3 x 10"" for AGS-791 ||. For purely 
left handed operators, the latter experimental limit places the bound A > 108 TeV. 
For scalar operators as in Eq. |3.1| , we find from Eq. A > 420 TeV. These results 
can be interpreted as bounds on the mass of new particles in different models. For 
example, we identify A y^rriH', the mass of a AG = horizontal gauge boson as 
discussed in Ref. [|^, to obtain rriH > 77 TeV. The bound on the family replication 
model of Ref. is placed by the scalar operators. Using Eq. ^]4| and g^l 
4v^(7^^ sin^ a/Ap we find Aq > 10^ TeV when we use the parameters sin a = 0.04 
and Qhc = 1- Similarly, we can use Eq. |3]^ to compare with the scalar leptoquark 
operators of Ref. |^ to obtain uih > 15 TeV, and mp > 2.8 TeV for pseudoscalar 
leptoquarks. In this case the bounds on the particle mass are not as strong as implied 



by Eq. because the couplings Cs±p in these models are suppressed by small mass 
ratios mf/mw- 

One can also use the form of the new operators Eq. p.l| to compare the reach of 
the rare kaon decay experiments with that of rare B decays. For example, taking a 
purely {V — A) (V — A) operator one finds 0: 

In the last step we have assumed that the B and K coupling constants are of the 
same order. This indicates that rare kaon decays are more likely to be sensitive to 
this kind of new physics than rare B decays. Of course, the situation may be different 
for other operators. 

These decays are also allowed in minimal extensions of the Standard Model in 
which the neutrinos are given a mass. The decays then proceed via box diagrams as 
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in Fig. ^ The decay rate for fi^e"^ is proportional to the product of mixing 

angles between the /x, e and the heavy neutral lepton A^: lUNeU^^l"^. Using the result 
B{fi 67) < 4.9 X 10-11 |§, the authors of Ref. |3 find |f/7Vef/^^P < 7 x 10^^ 
for m^v > 45 GeV. From this they conclude that B{Kl /^^e^) is at most 10"!^ 
in this type of models. Marciano has pointed out that there is a better bound 
If^TVe^TV/iP < 10^® coming from /i — > e conversion in the field of a heavy nucleus. 
With this bound one finds B{Ki /^^e^) to be at most a few times 10"!^. By using 
the additional theoretical prejudice that the mixing angles should be proportional to 
mass ratios of the form me/m7v, the authors of Ref. |^9|, |50| find the much stronger 
limit B{Kl — ^ /^^e^) < 10^^^. At this level the decay is completely unobservable. 
If, in addition, one considers left-right symmetric models, the authors of Ref. |4^, Q 
find that the rate B{Kl ~^ /^^e^) can be as large as 10"^^, although this happens 
only in a small corner of parameter space. 

There is an ultimate background to the decays yU^e^. It is given by the 

standard model process Kl fi'^e'^u^Ue. To leading order in xPT this process is 
dominated by Kl n^e^Ue followed by vr^ fJ'^^'^fi- The branching ratio for this 
chain can be estimated using the narrow width approximation to be a huge 38%. 
However, it is easy to see that the maximum invariant mass of the lepton pair is 
m^e < 489 MeV. It is therefore possible to remove this background with a cut on 
the lepton pair invariant mass. Going beyond the narrow width approximation, and 
including next to leading order terms in xPT can yield a lepton pair invariant mass 
larger than 489 MeV. However, after the 489 MeV cut is imposed, the rate is reduced 
to the 10-23 level \M. 
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3.2 K+ 7r+/i±e^ 

In this case we compare the rate to that of -k^ ^^v. From Eqs. |2.16| , p.l7| we 

can see that this mode is only sensitive to the vector quark current or scalar density. 
In general, this mode is thus probing different operators than the previous one. For 
V — A operators we find: 



(3.6) 



Fig. ID c shows the data from the most recent search for rc^fi'^e^ |^2|. The 

lack of candidates within the search region allows the 90% c.l. upper bound B{K^ 
TT^fi~^e') < 2.1 X 10~^° to be set. The corresponding bound on the state with opposite 
lepton charges, B{K^ n^fi^e'^) < 6.9 x 10^^ ||53[, was set in an earlier experiment. 
Using these experimental bounds we find A > 76 TeV. Once again, we can use this 
result to place bounds on masses of new particles within specific models. 

Given the large number of parameters in extensions of the standard model, it 
is very difficult to make definite predictions for any of the rare decays. Most of 
these models, however, will give rise to both K ^ fie and K vr/ie, as well as 
additional contributions to processes that do occur in the minimal standard model 
such as K ^ nuV, Kl and AM{Kl — Ks). Discussions of how to use the 

experimental limits on all these processes to constrain the parameters in the models 
can be found in Ref. [Q. 
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4 Short distance dominated processes 

These are modes that can be computed rehably from the standard model, because 
they are not affected significantly by the non-perturbative aspects of the strong inter- 
actions. Within the standard model they are sensitive to the CKM parameters p and 
?7, and to the top quark mass. Physics beyond the standard model, as in the models 
that gave rise to the modes studied in the previous section, can also affect these de- 
cay modes. We will concentrate on the standard model, but, as we will see, current 
experiments have not yet reached the sensitivity required to observe the standard 
model rates. This means that for the time being, these modes constitute a window 
to possible new physics. 



4.1 TT^UTJ 



In the standard model this process is mediated by the electroweak penguin and box 
diagrams depicted schematically in Fig. ^ The top-quark contribution is calculated 
to be 0: 

G p 



M = A'X%1- p-iri) 



a 



y/2 277 sin^ Ow 



X{xt) < TT+\s^^d\K+ > T7^^{1 + 75)z/ (4.1) 



where Xt = mf/ and 



X{xt) = ^ ( ^ + -^^^Inx,) (4.2) 
^ ' 8\xt-l {xt - 1)2 7 ^ ' 

Ref. [0 gives the approximate result X{xt) ~ 0.650x°'^^. The hadronic matrix 
element can be related by isospin to that occurring in —* Tc^e^u (equivalently 



using Eq. ^^.16] ), and this allows one to write: 



B{K+^Ti^e+p) VTTsin^^ 



a 



;i-p)'+^' 



X{xt) 



(4.3) 
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for each neutrino flavor. The contribution of the charm quark intermediate state 
has also been computed. In this case, however, additional kinematical factors for 
the lepton masses, as well as QCD corrections do not allow us to write a simple 
expression. The authors of Ref. [|TI| provide us with an approximate expression for 
the total rate with three generations of neutrinos: 

B{K+ 7r+z/I7) ^ 1.97 X A^xl'^^'lrf + (p^ - pf] (4.4) 



The parameter po measures the charm quark contribution and is given in Ref. |TT]|, a 
typical value for it being 1.5. This same reference finds B{K^ -k^vV) to be in the 
range (0.5 — 8.0) x 10^^*^ when all uncertainties are included. 

One can think of long distance contributions to this process that occur via /i pole 
diagrams as in Fig. ^, but these have been estimated to be at the level of 10~^^, much 
smaller than the short distance contributions 0. 

It has been argued that minimal extensions of the standard model are unlikely to 
affect this process significantly, so that its main interest remains the constraining the 
CKM parameters |^. However, even in Ref. ||55||, several less minimal extensions were 



discussed in which the current experimental upper bound could easily be saturated. 
In addition experiments sensitive to K'^ n^uV can also see n^X^ where 

X° is a new light weakly-interacting particle. Candidates for such an X° include 
axions [Q, familons [^, and hyperphotons [Q. Thus the possibility of novelty in 
the topology ^ 7r++ "nothing" should not be forgotten. 

Fig. 1^, shows recent data from AGS- 787. There were separate analyses for the 
kinematic regions with vr"^ momenta above and below that corresponding to the back- 
ground n^Ti^ process. No candidates were found within the search regions. 
The 90% c.l. upper limits extracted were B{K+ -n+uV) < 7.5 x 10~^ and 1.7 x 10"^ 
respectively p9|, |60| . Assuming a vector momentum spectrum for the tt"*", these can be 
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combined into an overall upper limit B{K^ — > ti^vV) < 5.2 x 10^^. The correspond- 
ing limit on 7r+X° for a massless, long-lived X° is 1.7 x 10^^. An upgrade of 
this experiment expects to reach a sensitivity of 10~^°, at which point one can begin 
to extract useful constraints on p and r). 



4.2 Kl^ TfOz/F 

The neutral version of the previous mode is completely dominated by the top-quark 
intermediate state. The result for three neutrino species can be written as ||TT|: 



B{Kl TT^uV) 3TK^f a ^2 



X{xt) 



2 



(4.5) 



B{K+ -K^e+u) 2 tk+ Vtt sin^ Ow ^ 
The expected rate in the standard model is of order 10^^°, and this decay would 



directly measure rj The potential long distance contributions to this decay are 
expected to be very small. Recently, the QCD corrections to the top-quark contribu- 
tion to this decay (and to the charged mode as well) have been computed in Ref. [|62 |. 



This calculation reduces the uncertainty in the rates due to the dependence of the top 
mass on the renormalization scale from (9(25%) to (9(3%). This makes this mode a 
particularly clean one to measure rj within the standard model. In principle this mea- 
surement can also be combined with that of the charged mode to reduce uncertainties 
due to rrit and to A. 



At present the best limit comes from FNAL-731 it is B{Kl — > vr'^z/I/) < 
2.2 X 10~^. In this experiment the decay was sought in the Dalitz mode where the vr" 
decays to e+e~7. Using this technique, FNAL-799 expects to reach sensitivity in the 
10~* region. The KAMI conceptual design report |Q claims an eventual sensitivity 
of a few xl0~^2 for this decay at the FNAL Main Injector. There is also a letter of 
intent for an experiment [^] at KEK with an eventual sensitivity of < lO^^^/event. 
Other versions of this experiment are under consideration at a number of institutions. 
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5 Long distance dominated radiative decays 

In this section we consider radiative decays and some of the Dahtz pair modes as- 
sociated with them. These decays are sensitive to the non-perturbative aspects of 
the strong interactions, and thus they cannot be predicted rehably from the standard 
model at present. We resort to the framework of xPT to study these modes. As we 
will see, the current state of the art, 0{p'^) calculations, is not always sufficient for 
an adequate description of these modes. These modes are interesting in their own 
right because they yield information on the long distance strong interactions. In some 
cases, a detailed understanding of these modes is also necessary in the analyses of 
other modes that look for CP violation, CKM angles, or new physics. 

5.1 Ks 77 

The processes Ks — > 77 and — > 7r''77 share the remarkable feature of being 
independent of the couplings in Eqs. |2.12| , |2.13| at 0{p'^). They are therefore among 



the cleanest predictions of xPT. They can be studied by considering the diagrams 
shown in Fig. ^. A considerable simplification in the calculation is obtained when 
one uses the "diagonal" basis of Ecker et. ai, Ref. Defining z = q^/^'x-, 

rl = ml/mj^, we write the amplitude for 7(q'i)7*(q'2) as restricted by gauge 

invariance: 

Mu = e^gi) (^(^ - l)9,u + Wi.)&(0, z). (5.1) 
The Oi^p"^) result is given by p8| : 

6(0, z) = -G82v^A(l - rl)H{z). (5.2) 

71 

The function H(z) can be written as: 



(5.3) 
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the function F{x) is given by: 



F{x) 

and G{x) is given by: 

G{x)- 



l-l(arcsinf)^ 



1 + iflogi^^ 



^-ifarcsin^y 



X < 4 



X > A ] 



X < A 



(5.4) 



(5.5) 



The situation in K^. 



77 looks very good, with 2; = in Eq. [5.2| , one finds a 
branching ratio of 2.0 x 10"^ [§3 and the NA31 measurement is (2.4±1.2) x 10"^ 
However, we note the very large error in the data. Given the fact that the theoretical 
prediction is rather clean, it is very important to reduce the experimental error. It is 
also possible to study this prediction in the related Dahtz pair decays ^"^£"7, 



for which one finds 28 



— = — r(ir°^77)-(i 

az 67r z 



Hiz) 



HiO) 



'l + 2Tl/z)Jl-Arj/z 



Ignoring CP violation the predicted branching ratios are then: 

T{Ks ^ £+£-7) _ [ 0-016 ^ = e 

T{Ks ^ 77) \ 3.75 X 10-4 ^ = ^1 

5.2 Kl 7r%7 



(5.6) 



(5.7) 



The amplitude for Kl{p) -> vr°(p')7(gi, £1)7(52, £2); M = ei(gi) €2 (52)^^1,, is re- 
stricted by gauge invariance and CP conservation to be of the form: 



M^^ = A{z,u)(^-m%^g^^ + q2f,qiu^ + 



-mj^XiX2g^,u 



where Xi = p ■ qi/mj^, z = 2qi ■ q2l'm\, and u = Xi— X2- The lowest order amplitude 
Kl — > 7r''77, C^(p^), obtained from chiral loops gives rise only to the form factor 
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i/)|^6[. To this lowest order result we can add the next order terms, 0{E ), as 
they appear in some models. They can be parameterized in terms of a single constant 
ay, and the combined result is: 



A{z,u) = -Gg 
vr 



r£/ V z 



F[-][l-^]+Fiz)[^-^-l]+ay[3 + rl-z 



a 

B{z,u) = -Aay-Gs . (5.9) 
vr 



The function F{x) is given in Eq. 

The constant ay has been calculated in several models. The simplest ones are 
those that consider only pole diagrams for the terms, such as those of Fig. |^. One 
should also consider direct weak counterterms as depicted schematically in Fig. |^. A 
model to compute these direct counterterms is the "weak deformation model" of [p6| . 
For this mode, the model predicts the direct weak counterterm contribution to ay 
to be twice as large as that from the pole terms and to have the opposite sign. The 
net effect is thus to change the sign of the constant ay calculated from pole diagrams 

0.32 VMD (poles); -0.32 (poles + WDM) [36] 
ay = { (5.10) 
0.22 Q.M. (poles); -0.22 (poles + WDM) 



We compare these results with experiment ||70|, |71| in Fig. |[ In Fig. || a, the 
various predictions have been mutually normalized. It is apparent in Fig. ^ b, that 
the shape of the data of Ref. ||70[ is well fit by the O(p^) prediction (this is also true 
of the data of Ref. ||TT|). The shape of the distribution is affected very little by the 



small corrections. The corrections (both in the quark model and in the vector 
dominance model) tend to increase the number of events to be expected for the lower 
values of the photon pair invariant mass. However, this is a very small effect, and it 
is difficult to check given the currently available statistics. 

The branching fraction with a cut on the invariant mass of the photon pair. 
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M^^ > 280 MeV, is 0, 0: 

B{Kl ^ vr%7) = <! 



f (0.57l°:J^) X 10- 



(5.11) 



(1.7 ± 0.3) X 10"*^ NASI, 

^ (1.86 ± 0.6) X 10-*^ FNAL-731. 
The central value for the theory number corresponds to the 0{p^) result. The error 

represents the variation obtained by including the 0{p^) terms with ay from Eq. ^.10 



The predictions are significantly smaller than the data. Other models (outside the 
realm of xPT), have predicted larger branching ratios |72|; however, they predict a 77 
spectrum that strongly disagrees with experiment, whereas xPT appears to reproduce 
the data very well (see Fig. H). These other models include vector meson exchange 
diagrams such as those of Fig. ^. This has created the incorrect impression in the 
literature that it is the inclusion of vector mesons that is responsible for the large 
rate. That this is not so can be seen, for example, in the model of Ref. Eq. [5.10 , 
that includes the vector mesons as well. As noted in Ref. p^, what leads to large 



rates in the models of Ref. |72] is the very specific form in which they include the rj' 
pole. This pole, unlike the 7r° and r] poles of Fig. |^, cannot be treated unambiguously 
at present. It enters at the same order asrj — r]' mixing and other higher order SU (3)v 



breaking effects. The authors of Ref. |]72| include only some of these effects, and the 
rate is very sensitive to the precise way in which this is done. The very large model 
dependence introduced by these terms can be appreciated in the mode Kl — * 77 



where a similar situation occurs 73 



From the shape of the spectrum, NA31 has derived the bound [[70[] : 



0.32 <av < 0.19 



(5.12) 



at 90% confidence level. One must notice, however, that with this range for ay, the 
branching ratio (with M^j > 280 MeV) is not reproduced: 



0.53 X lO"*^ < B{Kl n^ii) < 0.68 x 10"^ 



(5.13) 
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Conversely, a fit to the branching ratio yields ay ~ —2.0, which does not reproduce 
the shape of the spectrum. The shape of the spectrum appears to indicate that the 
xPT arguments are correct, and that models with very large D-wave amplitudes are 
ruled out. The discrepancy in the branching ratio tells us that this mode will require 
more theoretical and experimental efforts in the years to come. 

The contribution of the AI = 3/2 amplitudes has been recently estimated in 
Ref. [j74|. These authors found that it introduced a small (9(10%) correction that 
tends to make the branching ratio smaller. The amplitude computed in Ref. 



contributes only to the A{z, v) form factor, so the spectrum still has the shape given 
by lowest order xPT. These same authors have studied possible effects of higher 
order terms by using the experimental amplitude for the — > Svr vertex. They 
find an enhancement of about 26% in the branching ratio, that comes mostly from 
contributions to A{z^ v). Although this does go in the right direction, it is too small 
to account for the discrepancy with experiment. We should also remember that the 
naive estimate for the size of corrections to a amplitude is about 25% as well. 



5.3 K 



This process is similar to the neutral mode in that the one-loop amplitude in xPT 
is finite. Unlike the neutral mode, however, this mode does receive contributions 
at order 0{E'^) from local counterterms. An additional feature is the possibility of 
looking for CP violation by comparing the two charged modes. The amplitude for 
this process is restricted by gauge invariance to be of the form: 



M 



ei e 



(5.14) 



The form factor A{z^ v) to 0{p^) is given by p8[] 

z 



A{z,v) 



a 1 

TT 2Z 



z ]F 



+ 1 



ri]F(z) 



cz 
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c = 327r'(^4(L9 + Lio)-^(wi + 2«;2 + 2w4)) (5.15) 

The constant c is a scale independent combination of couplings in the effective La- 
grangian, reflecting the fact that the loop contributions to this process are flnite. Its 
real part, Rec, has been estimated in several models to be: 0.0 |3^; 0.35lo;3o 
and -4.0 p7|. 



The form factor C(z, z/) occurs via the pole diagrams of Fig. ^. The transition 
71^71^ occurs on-shell via the (27l, 1r) operators only. The 7r+7r° decay, 

followed by tt^ '-ff gives a very large contribution to 7r+77, however, it 

can be subtracted by implementing a cut in the invariant mass of the photon pair 
to exclude the region near m^o. The off-shell transition is mediated by the octet 
operator in Eq. and results in [p8| : 



C{z, u) = -G, 

71 



z-rl z-|(2 + r2) 



The two form factors do not interfere, and one obtains: 



(5.16) 



dT 



,5 



:z'xHl,z,rl)(\A{z)\' + \C{z)\'^ (5.17) 



dz 10247r3 

In this case, xPT predicts a correlation between the rate and the spectrum. The rate 
is given by the expression: 

B{K+ ^ 7r+77) = (5.26 + 1.64c + 0.320^ + 0.49) x 10"^ (5.18) 

and the current experimental limit is: 



1.0 X 10"^ Phase space spectrum 
B{K+ ^ 7r+77) < { (5.19) 

1.5 X 10^^ xPT spectrum 

The rate for K~ 7r~77 can be obtained from Eq. ^.15| by replacing Gg and c with 

their complex conjugates. An imaginary part of c would interfere with the absorptive 

part of the amplitude (due to a real tt^tt^ intermediate state) generating a CP odd 
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rate asymmetry [p8|: 



r{K+ 7r+77) - r{K~ 7r"77) = 1.5 x lO^^^Imc GeV 



(5.20) 



To estimate Imc, the authors of Ref. point out that at the quark level, the CP 



phase appears via the electromagnetic penguin operator |76|: 



£ = -r^|KrfK.|C7(/x')as7^(l + 75)dir^ 



(5.21) 



This operator transforms as an octet under SU{3)v- By requiring its chiral realization 
to transform in the same way, the authors of Ref. conclude that: 



(5.22) 



Arguing that this is the only contribution to the phase of c they find: 

32 

|Imc| = — TT^IImwil ^ 0.00577 (5-23) 

Noticing that the rate for 7r+77 given by Eq. ^.1^ has a minimum for c = 

—2.56, they conclude that: 

T{K^ — *• 7r"'"77) — T{K~ — *• 7r~77) 



< 0.002r7 



(5.24) 



T{K+ 7r+77) + T{K- 71-77) 
However, the authors of Ref. have pointed out that the electromagnetic penguin 
operator also contributes a phase to W4. They find that this contribution essentially 
cancels the one from wi so that Imc = and there is no charge asymmetry in the 
standard model. Given the potentially large asymmetry, Eq. p.24| , it is important to 
resolve this issue. 



5.4 Kl '-f-f 

This process is not a very rare decay, but its calculation in xPT illustrates some of 
the problems that one encounters in other weak decays. To order 0{p'^) this process 



-28- FERMILAB-Pub-93/004-T 



occurs via the 7r°, t] poles of Fig. [1^. Using Eqs. ^]8| and 2]9| one finds in the hmit of 
CP conservation that M{Kl l{.<li)l{.<i2)) is: 

M = + l!4z4),_,.,J,;4 (5.26) 

At this order in xPT, however, we are instructed to use the Gell-Mann-Okubo mass 
relation Eq. p.6|, so the amplitude vanishes. The first non-zero contributions result 
from higher order S'f/(3) breaking. We can parameterize these higher order effects 
with a constant Fi, such that: 

M = -2G8AFie^,,.gfg2^e?e^ (5.26) 

TT 

to find: 

r (7.3 X 10-4)|Fi|2 
B{Kl 77) = <^ (5.27) 
[ (5.70 ± 0.27) X 10"^ PDB average [40] 

A fit to the experimental result yields a value ~ 0.88. A detailed analysis that 



tried to predict the value of Fi was carried out in Ref. [^. It included SU (3) breaking 
effects in the vertices as well as the decay constants; and also the effects of the SU (3) 
singlet and of rj — t]' mixing. Ref. |j73| found that it was possible to accommodate 



the experimental result, but that it was not possible to predict it with any certainty. 
Recently, the authors of Ref. |]7^ have criticized the treatment of the SU{3) singlet 
in Ref. |7^, however, this does not change the conclusion. 

It is also of interest to study this amplitude for off-shell photons, since this will 
give the dominant contribution to the decays — > and — > . In 

the limit of CP conservation, gauge invariance requires the amplitude for —>■ 'J*'J* 
to be of the form: 

M = MiKL ^ l{qMq2))Ciql, ql) (5.28) 

where we have normalized it so that C(0,0) = 1. Bose symmetry requires C{ql,q2) 
to be symmetric under qi ^ q2, so we write: 

Ciqlq!) = l + a(^^±^)+.-- (5.29) 



m 



K 
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5.5 Kl 77 Dalitz decays 



Using the previous result, and defining = me/rriK, z = mjjmj^, we can now study 
the process Kl — > The differential decay rate is: 



dr 

dz 



Y(Kl 77) 



2a 
3^ 



CiO,z] 



1 



4ri 



1 



(5.30) 



and the total rates are: 



r{KL 77) 



r 0.016(1 + 0.05a) 

0.016 ±0.001 

4.09 X 10"^(1 + 0.64(t) 
^ (6.81 ±0.64) X 10-^ 



AGS-845, NA31 [78, 



(5.31) 



FNAL-799 ISO 



It has become usual in the literature to parameterize the form factor following a 



model of Bergstrom et. al. fSlj], in terms of a constant ax- This model, with ax = 
corresponds to the vector dominance model of Quigg and Jackson |8^. Expanding 
the pole model parameterization, we identify a = 0.418 — 1.29q;a'- In terms of this 
pole model, AGS-845 [0, |7| found ax = -0.28 ± 0.083 by fitting the m« spectrum 
with QED radiative corrections, and ax = —0.18 ± 0.077 without QED radiative 
corrections. This is consistent with the NA31 result, ax = —0.28 ± 0.13. 

The branching ratio T{Kl e~^e~'~f)/r{KL — > 77) is not very sensitive to the 
form factor. With the form Eq. |5.29| , we find a variation from 0.0162 to 0.0165 when 



we take ax from to —0.28. With the pole model of Ref. the branching ratio 
varies from 0.0163 to 0.0167 for this same range of ax- 

The branching ratio T{Kl ^'^ ^^^) /T{Ki 77), is more sensitive to the form 
factor. With the same range for ax as before, it varies from (5.18 to 6.13) x 10""^ using 
Eq. |5.29| ; and from (5.58 to 7.28) x 10~^ with the model of Ref. jST]]. The preliminary 
result from FNAL-799 cited in Eq. |]31, can be used to extract from Ref. [p]|: 



ax 



-0.21 ±0.11 



(5.32) 
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Which is consistent with the value of ax found by the — > e"'"e~7 experiments, 
and inconsistent with = at the two a leveL Prehminary reports from FNAL- 
799 indicate that the m« spectrum for this decay is adequately fit by an unmodified 



Kroll-Wada form 33 



Finally, there have been several recent results on the double Dalitz decay, 



e^e e^e |83, 85, 86, 87|, which are summarized in Table H 



Table 3: Experimental Results on 



e ' e e ' e 



Experiment 


Result 


Comments 


NA31 


(4 ±3) X 10-*^ [34] 


the first 2 events 


AGS-845 


(3.04 ± 1.24 ±0.26) x 10"*^ |85| 


6 events 


KEK-137 


(6 ±2 ± 1) X 10-s [36 


partial reconstruction 


FNAL-799 


(4.17 ±0.83) X 10-s [0 


28 events, preliminary result 



These results should be compared with the theoretical expectation of 3.6 x 10~^ 
Larger samples will allow form factor effects to be studied. FNAL-799 also 
expects to observe the closely related process /i+yU^e+e^ which is predicted 

at 0.8 X 10-9 [H. Reference H also predicts B{Kl /i+/i-/i+/i-) = 5.4 x lO'^^^ 
putting it beyond reach for the moment. The predictions for these modes that one 
finds in the literature do not always agree [^ 



5.6 Direct emission K 



In the limit of CP conservation, this process does not occur at order 0(]5^) [^, pO [. 
The lowest order weak Lagrangian Eq. ^]9| contributes to the process 7r+7r^7, 
and thus it contributes to Kl n^n^^ an indirect CP violating term. This term has 



a characteristic bremsstrahlung spectrum and has been observed [0, |9^ . Of greater 
interest to us is the so called direct emission term. From a theoretical point of view. 
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we will define it as the amplitude that starts at 0{p'^). Experimentally, it is observed 
by subtracting the bremsstrahlung portion from the full amplitude. Gauge invariance 
requires the amplitude for Ki{k) — > Tf^ {p'^)tt~ {p~)'~f{q) to be of the form: 



M 



( 

rriK 



^e{z, v) f z/(j9+ +P )^,- z{p'' -p + -^iM{z, i')e^uapk''p'^°'p 
V / mi- 



(5.33) 

where z = 2k-q/m\, = 2k-p^/m'j^, and u = —x^. The "magnetic" form factor 
C,M receives contributions from 7r°, rj poles as in Fig. |ll|, and vertices from Eq. p.8| . At 
order 0{p^) we find: 



Cm 



m 



K 



m\ — ml Stt^/^ V 3 m\ — nii 



(5.34) 



However, when the Gell-Mann-Okubo relation (Eq. |2.6|) is applied, one gets the same 
cancellation that occurred in Kl 77. We can parameterize the SU{'i)v violating 
terms that make the amplitude non-zero, as well as terms of order 0{p^) with a naive 
expansion of the form factors. We follow Ref. and include the strong rescattering 
phases 5 j for the spin J, isospin / tt — tt scattering at a center of mass energy squared 

{p+ + P-Y = f^^i 

^M{z,iy) 



4(1 



z). With all this we write: 

v+- 



K 



mz 



+ 



m 



K 



l + aM;^)e^(^i~^o°)+,^,,^e^(4'-^E!) 



(5.35) 



We have included in the inner bremsstrahlung contribution which can interfere 
with Fe in the rate. We have not included possible 0{p^) AI = 3/2 terms. The 
terms Fe, and qm sue CP violating. The differential decay rate is given by: 

dT arriK 



dzdv 1287r 



z){z' 



Arlz^ 



(5.36) 



A detailed analysis of the Dalitz plot for this decay should allow a determination 
of the parameters in Eq. p.35| . The simplest thing we can do is to set all the new 
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constants in Eq. p.35| to zero except Fm- In that case we find for the direct emission 
(with E^ > 20 MeV): 



B{Kl TT+TT j)dE = < 



' (3.35 X 10-5)|Fa/|2 

(2.89 ± 0.28) X 10-^ BNL [91] (5.37) 
^ (3.19 ± 0.16) X 10-^ FNAL-731 [||] 
which results in Fm ^ 0.95. One can try to predict Fm in a model that includes 
SU (3) breaking effects, similar to the one used by Ref. [^] for — * 77. One finds 
that it is possible to accommodate the result [^ , but not to predict it with certainty. 
If we set Fm to zero and keep Fe in Eq. |5.35| , we find that this direct emission rate 
is also consistent with an interference of the bremsstrahlung amplitude and a CP 
violating term Fe ~ 1.4. Keeping only Fm and Fe, Ref. [^ found a best fit to the 
spectrum with an admixture of CP violating direct decay Fe ~ 0.26Fm, lying < 2cr 
from the no- interference fit. However, naive dimensional analysis suggests that CP 
violation in the standard model is much smaller than this [|93[ . The fits of Ref. [^ 
found no evidence for the quadrupole term qe at the 25% level. 

Assuming CP conservation in the direct emission, a fit to the distribution dV/dz 

provides information on the slope cta/ : 

f -0.91 FNAL-731 
<yM = (5.38) 
I -0.89 BNL 



One can resort to models to predict the 0{p ) form factors like ctm [p3| , [9^ , |95| . It 
appears that a vector meson dominance model does not reproduce the spectrum shape 
unless one includes SU{?>)y breaking effects. A more careful study is needed. 

By observing the interference between Ks and Kl decays into 7r+7r~7, FNAL-731 
9^ has measured: 



= ^}^' "^"^^ = (2.15 ± 0.26 ± 0.20) X 10"^ (5.39) 
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with a phase = (72 ± 23 ± 17)°. Since |?7+_| = (2.268 ± 0.023) x lO^^ gg, the 
FNAL-731 resuh is consistent with no CP violation beyond that present in K ^ titi 
decays. 

Additional information can be obtained by studying the Dalitz decay Kl — > 
7T~^ TT~ e" . Separating the I.B and D.E. contributions, and using for the D.E. a 
constant coupling Fm that fits the {Kl t^^t^^ i)de rate, Ref. ||9^ calculates 



B{Kl —>■ TT+TT^e+e^) = (1.3 + 1.8) x 10^''. This decay has been recently observed by 
FNAL-799, but a branching ratio has not yet been reported. The authors of Ref. |^ 
have also studied a CP odd correlation between the vr+vr" and e~^e~ planes and pre- 
dict a 14% asymmetry. Unfortunately this asymmetry is mostly due to indirect CP 
violation through the parameter e, making an extraction of direct CP violation very 
difficult. 
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6 Decays into charged lepton pairs 

These modes generally receive both long and short distance contributions that are of 
comparable size. Some of them are dominated by the long distance component, but 
exhibit interesting interference effects with the smaller short distance contribution. 
Their primary interest is the study of the short distance parameters p and 77, but 
to do so one needs some understanding of the long distance contributions. In this 
respect they also allow one to test the ideas of xPT. The observables that probe the 
short distance physics are also sensitive to new interactions, however, we will restrict 
ourselves to a discussion of the standard model. 



6.1 Short Distance Kl 

The short distance contribution to these processes comes from the box and electro- 
weak penguin diagrams of Fig. |12|. To compute the amplitude for this process, one 
relates the matrix element < 0|s7^(l + 'y^)d\KL > to that occurring in £+z/ 
(see Eq. p. 161 ). The top-quark contribution is easily computed to be |I0|] : 



B{KL-^i+i-)sD TK.f a 



where the function Y{xt) is given by: 



)W(l-p)^ 



Y{xt) 



2 



(6.1) 



An approximate expression is provided by Ref. [P] Y{xt) ~ O.Sldxf^^. The contribu- 



tion of the charm-quark, with perturbative QCD corrections can be found in Ref. [11 



These authors provide us with the approximate expression for the complete result: 
B{Kl ^ fi^fi-)sD = 1.7 X 10-i°x,i-^^A^(po - pf (6.3) 

where deviations of po from 1 measure the charm-quark contribution with QCD cor- 
rections. For typical values of all the parameters involved, po ~ 1.27. 
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An interesting feature of this decay is the longitudinal polarization of the final 
lepton, a CP violating observable [ P^ : 

Nt - Nr 

The latest estimate of this quantity within the standard model is about 2 x 10^^ 



This number, however, is directly proportional to e [Q. Since this is probably too 
small to be measured, Pl is a very good place to look for direct CP violation outside 
the standard model. 



6.2 Long distance t^t 



The long distance contribution to these decays is expected to be dominated by the 
two photon intermediate state as in Fig. |T3|. It is straightforward to compute the 
absorptive part of the amplitude by using the experimental rate for Kl —>■ 77. The 
result is (r| = mj/m'j^): 



BiKr 



1 2 2I 

-a r.^ 



In 



l + f5 



B{Kl ^ 77) 



1 - 4r? 



This can be compared to the latest measurements: 

' (6.8 ±0.3) X 10^9 



B{Ki 



B{Ki 



e e 



absorptive 

(7.9 ± 0.7) X 10-9 KEK 137 [100] 

, (6.86 ± 0.37) X 10-9 AGS-791 ^U^ 
(3.0 ± 0.1) X 10-^2 absorptive 

< 4.7 X 10"^^ 



AGS-791 42 



(6.5) 



(6.6) 



A calculation of the dispersive part of the long distance effects is not possible at 
present within xPT. There are two problems: the on-shell 77 vertex cannot 

be computed reliably as already explained; and the divergence of the loop diagram 
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implies that there are counterterms for this process. These counterterms are similar 
to those occurring in the decay 7r° 



e ' e 



102| |, but they are not the same and have 



not been determined. However, we can resort to models to estimate the dispersive 
contribution. A vector dominance model gives results that can be written compactly 



in the limit m^/mp — ||103 



2^2 r 



r(po 



X 



2a m 



nil 



+ In 



mi ) 



2n 



TT 



, ^ + irin!^V--ln— 

12 An TT V m£ / n mi 



(6.7) 



where my is the rho mass. Although this is a convenient expression for decays into 
electrons, it is not accurate by about a factor of two for muons. One also finds in the 
literature the model of Ref. in which one allows an additional form-factor for 
the Kl — ^ 77 vertex itself. This is modelled by a K* p transition in terms of the 



parameter ax of section |5.5| . The result can be written as [^T 
B{Kl fi^p )dLp 



-4.7x 10-^° ( 1.3 + 4.9a/^ 



(6.8) 



and the ax = case corresponds to the vector dominance model of Quigg and 
Jackson. The contributions from a negative ax tend to cancel the first term in the 
above result, and this cancellation is almost complete for ax = —0.28. We regard this 
cancellation as accidental, and prefer the result with ax = as a more conservative 
estimate for the long distance dispersive rate. 

In view of the large model dependence in estimating the long distance contribution 
to the dispersive part of it would seem very difficult to extract mean- 

ingful constraints on the short distance contribution. Nevertheless, recent attempts 



to do so can be found in Ref. [104 



From Eq. |6.7| we can also compute the two-photon contribution to B{Kl 
e~^e~)disp, we find 4.6 x 10^^^. The long distance contributions to Kl e^e~ are 
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about 15 times as large as one would expect by a naive scaling of — > ji^ ji^ with 
ml/rn^^. This factor can be traced back to the logarithms in Eqs. |6.7| . There 
is no analogous factor in the short distance contribution, which then has the naive 
scaling. This means that it is even harder to observe the short distance contribution 
in Kl e~^e~ than it is in Kl jjr 0]. 



6.3 K 



This is the most important long distance contribution to the CP conserving decays of 
the form K — > ni'^i^ . It can be computed to O(p^) in xPT in terms of a few unknown 
constants. The most general form allowed by electromagnetic gauge invariance for 
the amplitude A{K^'°{k) 7r=^'°(p)7*(e, g)) is 

eG, 



(6.9) 



where z = /m\. 

It can be seen immediately that gauge invariance requires the amplitude to have 
at least three external momenta (and there is one external photon), so that it can 
only start at 0{p'^). Eq. |6.9| also shows that the process K TC'-f with an on-shell 
photon is forbidden by gauge invariance. Computing the diagrams in Fig.p 



one 



finds the form factor C{z) to 0{p^) in xPT in terms of the couplings of section ^ and 
the function: 



- 1 ) atan,/^ - - + I 



z-2-y/z{z-4) 
z-2^^J7{I^ 



-, ~ e. 



z<A 
z > 4 



(6.10) 



The function (pT^^z) is obtained by replacing z 
the charged mode p7[]: 



zmj^ / m\ in Eq. |6.10| . One finds for 



(47r)^ <(/x) - wl{^,) + 3(4vr)^U;2^(/i) - ALl{^,) 
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+ In 



- (pAz) + (Pk{z, 



w+ + (f)^{z) + (Pk{z) 



and for the neutral mode: 



V2 
6 



(6.11) 



(6.12) 



6.4 K+ 



TT 



-i+t 



We can now study the processes K~^{k) 'K^{p)C^{k^)i {k ) that are dominated by 
the one photon intermediate state. The matrix element is given by {z = {k—pY/rn^^): 



47r 



C+{z)u{k-){f + i))v{k^ 



(6.13) 



From this, it is straightforward to compute the decay distribution p7| : 



.rj\2 



(6.14) 



dz 127r(47r)4' 

This is shown in Fig. |15|. Integrating Eq. |6.14| from 4r| to (1 — r-^Y and using 
Eqs. [6.101 , |6.11| we can get a prediction for the rates tt^ 
unknown constant 



in terms of the 



B{K+ ^ vr+e+e" 



n fi 



(3.15 - 21.1W+ + 36. Iw^) X 10" 
(3.93 - 32.7u;+ + 70.5w^) x 10" 



(6.15) 



One can then determine the constant from measurements of the rate, the spec- 
trum, or both. A measurement of w+ in vr+e^e^ has recently been re- 
ported [|105|| (AGS- 777). Eq. |6.16| and Fig. ^ give the result of a simultaneous fit of 
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the rate and spectrum shape: 

^ TT+e+e") = (2.99 ± 0.22) X 10"^ 

w+ = 0.89l°:?l (6.16) 



Eq. |6.15| , which has not been imposed on the fit, is shown as a parabola in Fig. |T6|. It 



passes within about 1.5a of the best fit values. As consistency checks, one can insert 
these values in turn into Eq. |6.15| yielding: 



B{K^ ^ TT+e+e-) = (1.30+JJ^) x IQ-' 

r 1.20 ±0.033 or 
w+ = (6.17) 
I -0.62 ±0.033 

This approach to 7r+e+e^ seems quite promising, but more data is desirable. 

A somewhat larger data set with better systematics should be forthcoming from AGS- 
851. In the longer term, a much larger sample is promised from AGS-865 which is 
presently under construction. 

Some model calculations of w+ have found: 0.7 |3^; 0.98iJ]:^^ |8[; and 1.9 g. 
If we take as a typical number = 0.89 we then predict B{K^ vr+yU+yU^) = 
3.07 X 10^*^. At present there is only an upper limit of < 2.3 x 10^^ on this decay 
from AGS-787 [ 106 |, but a data set with sensitivity < 10~^/ event is presently under 
analysis by this experiment. 

A CP violating imaginary part of the constant w+ would interfere with the ab- 
sorptive part generated by the two pion intermediate state (Eq. |6.10|) , giving rise to 
a CP odd rate asymmetry p8| : 



T{K+ vr+e+e-) - T{K- Ti-e+e 



0.01Imu;+ ^ 3 X 10"^//, (6.18) 



T{K+ TT+e+e-) ± T{K- vr-e+e" 

where in the last step we have used the estimate of Ref. [^, Imw+ ^ 0.003?]. This 
is too small to see in the near future. 
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We should comment on a possible parity violating asymmetry for these decays. It 
was pointed out in Ref. | p.07| ] that apart from the one-photon intermediate state, these 
processes have a short distance contribution from a Z intermediate state or from box 



diagrams as in Fig. |12[ These new operators generate a second possible amplitude in 
addition to Eq. |6.13 : 



(6.19) 



where we have used the lowest order xPT result /+ = 1, /_ = 0. The constant ^ 
contains the short distance factors, it is given in Ref. | fL08| |: 

Y{xt) 



i = -1.4 X 10"^ 



(6.20) 



2n sin^ 9w 

where Y{xt) is given in Eq. |6.2| . The first term in Eq. |6.20| corresponds to the charm- 
quark contribution with QCD corrections, and the second term to the top-quark 
contribution. If we denote by T^, T^, the rates for producing a left-handed (right- 
handed) fi'^ in K'^ 7[~^fi~^fi~, then the interference of the two amplitudes generates 



the parity violating observable ||108 

\Alr\ = 



R 



2.3ReC 



(6.21) 



A measurement of |Ai^ij| at the tenth of a percent level would provide valuable in- 



formation on the CKM parameter p | 108 , 109 |. Taking for example, rrit = 140 GeV, 
p = -0.51, Eq. yields | A^^] = 3.7 x lO'^. 

Finally, the authors of Ref. | |110| ] have proposed some T-odd observables that can 
be studied in this decay. Unfortunately, in order to extract information on CP vi- 
olation from this type of observable one must be able to reliably subtract unitarity 
effects. This usually involves a comparison of the two charge conjugated modes. In 



Ref. [110] it is found that the unitarity effects are small in asymmetries where the 
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polarization of both the /i"*" and the /i~ are measured. Unfortunately such measure- 
ments are extremely difficult, due to the absorption of stopped fi^ before they can 
decay. 



6.5 TT^^+r 



In the limit of CP conservation, the one-photon intermediate state contributes only 
to Kf{k) ^ ■K°{p)i+{k+)i-{k-). The decay distributions can be predicted in terms 
of the constant Ws with Eq. |6.14| simply using: 



C'^iz)=ws + 2Mz) 



(6.22) 



Ignoring CP violation this results in rates and spectra for the decays Ks —>■ ir^i^i 



For the rates one finds ||106|| : 



B{Ks 7r°e+e" 
B{Ks ^ 7r'fi+fi- 



(3.07 - 18.7ws + 28Awl) x 10 
(6.29 - 38.9^5 + 60.1wl) x 10 



-10 



-11 



(6.23) 



In general, xPT does not relate ws to w+ measured in the charged mode. Without 
additional input one must first measure in one of the decays and then use it to 
predict the others. The prediction one finds in the literature is based on an additional 
assumption. The authors of Ref ||28| demanded that the meson Lagrangian transform 
under SU{3)v as a pure octet in analogy with the quark electromagnetic penguin 
operator. That gave them the constraint W2 = 4:Lg which then allowed them to 
predict: 

(6.24) 



1, f mi^ 
ws = w+ + -In — 
6 \m 



K 



One must remember, however, that Eq. |6.24| is an assumption that goes beyond xPT 



and that it is not satisfied in some models |38|. Sample model calculations of ws 
yield: 0.3 |3; 0.98j:[}|^ |8[; and 1.4 [|36|. 
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6.6 Kl Tv^e+e- 

The decay Kl vr^e^e" is significantly more complicated than the others we have 
been discussing. It has at least three different contributions that could be of the same 
size. The most interesting one, of course, is the direct CP violation. It originates in 



the diagrams of Fig. O. The result has been computed in Refs. [11, 76, 111]. The 



full result has a complicated form, but Ref. [0| gives an approximate expression: 

B{Kl 7r°e+e") = 0.32 x 10'^'^7]'^A^I{xt) (6.25) 

with I{xt) ~ 0.73x1'^^. With present day bounds on all the parameters, Ref. finds 
that this contribution to the rate ranges from about lO"^^ to 2 x 10"^^ 

There is also an indirect CP violating contribution, that is, one that proceeds via 
the parameter e in the mass matrix. Its contribution cannot be computed directly at 
present, and its precise value will only be known after a measurement of Ks 7r°e"*"e~ 
is done. However, one can also use Eq. |6.22| to predict the indirect CP violation in 
the decay Kl irH+i--. 



B{KL^^^H^n^nd = \e\'^BiKs^7r'i^n 
B{Kl vr°e+e-)i„rf < 1.27 x IQ-^^ (6.26) 



The last result follows from using w+ = 0.89loj4, and Eq. |6.24| to obtain ws 



ai4- With this range for ws, the rate Eq. |6.26| varies by more than three orders 
of magnitude! In fact, for ws = 0.33, Eq. |6.23| is not sufficiently accurate to calculate 
the rate. Given this large sensitivity to ws-, and the fact that we rely on Eq. |6.24 , 
this result must be viewed with caution. 

Finally, we can use the result for Kl 71^77, Eq. |5]^, to estimate the rate for the 
CP conserving part of the Kl — > 7r°e+e^ amplitude. As usual ||112|| , we will simply 



give the contribution from the absorptive part of the two photon intermediate state. 
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depicted in Fig. |T8|. The contribution from A{z, v) is suppressed by rrie and can be 
neglected. Using Eq. |5.9| we find a simple result if B{z^ v) is constant or if it depends 
only on z. We find p^ : 



AaUKdv) ^ 7rOe+(A;')e-(A;)) = -^a^^-^p ■ {k - k')ui>v (6.27) 

After squaring and integrating over phase space, this gives a lower limit for the 
branching ratio from the CP conserving amplitude. With —0.32 < ay < 0.19, Eq. 



5J2| , Ref. quotes: 

Bcp{Kl 7r°e+e-) < 4.5 x lO'^^ (6.28) 

However, we must remember that this is only the absorptive part of the amplitude, 
so that the number is not really an upper bound. However, it is sufficiently smaller 
than the direct CP- violating component, that it will probably not impede efforts to 
extract the latter. 

The fact that the direct CP-violating contribution to this decay is comparable or 
greater than the competing contributions, has sparked a good deal of experimental in- 
terest, including a number of dedicated searches. The present situation is summarized 
in Table g 

Although this decay has a good kinematic signature and its all-electromagnetic 
final state can be exploited in the design of experiments, it has unfortunately been 
found to suffer from the complication of a very difficult background. This stems 
from the processes shown in Fig. |19|, i.e. radiative corrections to —>■ e~^e~'j, 



resulting in e^e 77. This was first observed by AGS-845 [|116|| in the course 



of their search for Kl — ^ vr^e+e [|115|| . The branching ratio was measured to be 



(6.6±3.2) X 10 ^. The potency of this process as a background to Ki 7r°e+e , which 



had not previously been appreciated, was explicated by Greenlee |117|. Although 



there is no particular enhancement near m^j = m^o, he found that the rate is sufficient 
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Table 4: Summary of — > Tr^e^e Experiments 



Experiment 


Kesult 


Status 


Comments 


1\T A Q 1 

iN Aoi 


< 4 X iU iio| 


nnisned. 




FNAL-731 


< 7.5 X 10^9 [ii4| 


finished 




AGS-845 


< 5.5 X 10-9 |115| 


finished 




FNAL-799 




analyzing 


aims for lO^^*^ — 10^^ sensitivity 


KEK-162 




running 


aims for 10~^° sensitivity 


FNAL-799II 




under construction 


aims for < 10~^° sensitivity 



to make the extraction of a — > 7r°e"*"e~ signal at the 10~^^ level problematic. For 
example, assuming an acceptance interval Am^^ = 5 MeV, with highly optimized 
kinematic cuts that accept half the —>■ -K^e^e' events, the background enters at 
an equivalent branching ratio of 10~^°. If the cuts are further tightened to include only 
10% of the signal events, the background is reduced by only about a factor 3. Further 
progress can be made if the resolution in m^^ can be improved beyond the already 
very optimistic assumption used by Greenlee. Another approach to coping with this 
background is to attempt to subtract it, capitalizing on its smooth dependence upon 
m^^. Of course this requires high statistics, which are not easy to come by at the 
required level of sensitivity. 

6.7 Kl TT^fl+fl- 

This decay mode, unlike the previous one, receives a substantial contribution from 
the two photon amplitude at 0{p'^) since the muon mass is not negligible. This results 
in the possibility of substantial interference between the CP conserving and violating 
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amplitudes. The amplitude can be written as: 



M = —ReGsu{k-)\imMz) - {It + :6)g{z)]v{k+) 



The CP violating form factor is 



g{z) = e[Rews + 20i^(z)] + ilmws 



(6.29) 



(6.30) 



whereas the CP conserving form factor is given by ||28 

1-/? 



h{z) 



a 



[z - ri)F 



z-l-ri)F(z) 



(6.31) 



where (3 = — Ar^/z. An interference between the two amplitudes generates a 
CP violating muon polarization. Taking for example Rew^ = 0.73 and Imws = 0.001 
Ref. [0 finds an average transverse muon polarization < ^ >= —0.37 and a branching 
ratio B{Ki n^^^fi^) = 6.3 x 10^^^. At present there is a preliminary limit 
B{Kl ttVV") < 1-7 X 10-*^ from FNAL-799 g. 
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7 The experiments 

In this section recent and current experiments are described and their results summa- 
rized. We exclude those results not directly related to the subject of this review (e.g. 
e'/e, rare vr decay). We also discuss briefly new experiments now under construction. 

Although these experiments span a range from to greater than 100 GeV/c in 
beam momentum, they have certain important features in common. In each case 
the source of the kaons is the interaction of protons from a synchrotron with a fixed 
target. An intense secondary beam is created and transmitted to a decay region 
viewed by a detector. In most cases the kaons in the beam are outnumbered by other 
species of particles. These, along with the kaons which don't decay, have to either be 
transmitted through insensitive regions of the detector or somehow absorbed without 
doing irreparable mischief. All but one of the experiments use magnetic spectrometers 
to measure the momenta of the charged decay products. The neutrals are measured 
in calorimeters of various types. Scintillator hodoscopes and particle identification 
devices such as atmospheric Cerenkov counters and muon filters provide triggering 
capability. 

The experiments all exploit common K decay modes for calibration and normal- 
ization. Very often the same modes used in this way are sources of the backgrounds 
that have to be confronted. These backgrounds must be fought both at the trigger 
and analysis level. Particle identification, timing, geometrical, and kinematic selec- 
tion reduce the large data samples collected to manageable size. Except in one case, 
the signal events are completely reconstructed. Typically the last stage of the analysis 
is a two-dimensional plot of effective mass of the final state particles versus a variable 
which reflects their direction. The signal is sought in the region of the K mass and 
small angle with respect to the beam (see e.g. Fig 
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All the experiments study at least one "tune up" process on which to demonstrate 
their ability to detect rare decays. This is a decay topologically similar to the primary 
object of the experiment, but somewhat more copious. These are often of considerable 
interest in themselves. For experiments seeking fie there is — ^ /W^/W ; for 

TT+yU+e^ there is — > Tr+e^e^; for ir^uu there is 71^77; for 

Kl TT^e'^e~ there is Kl — > 7r''77 or Kl — > 7766. 



7.1 AGS-845 



Fig. |2y is a plan view of the apparatus designed by a BNL-Yale collaboration to 
perform the world's first dedicated Kl TT^e~^e~ experiment, AGS-845. It was 
optimized to detect all-electromagnetic Kl decays (the lead filter and hodoscope at 
the rear were used to veto penetrating particles). Several million Kl (along with 
~ 3 X 10*^ neutrons) entered the 6-meter evacuated decay region during each 1-second 
AGS spill. A single-magnet drift chamber spectrometer measured momenta and a 
lead glass Cerenkov array detected 7s. The latter also served to measure the energy. 
Comparing this energy with the momentum distinguishes these particles from 
pions and muons. A 2-m long atmospheric hydrogen Cerenkov counter completed the 
particle identification. 

The expected potential backgrounds to Kl —>■ 7r'^e+e~ were Kl 27i^ or 3tt^ in 
which two of the vr'' undergo Dalitz decays, and accidental coincidences of 2 7s with 
Ke3 decays wherein the vr is mistaken for an electron. AGS-845 was able to eliminate 
all background and set a 90% c.l. limit of B{Kl 7r*^e+e^) < 5.5 x 10~^. This 
represents a large improvement in our knowledge of this process, but falls short by at 
least two orders of magnitude of the Standard Model prediction for this CP- violating 
decay. However, as discussed in subsection , in the process of setting this limit. 



AGS-845 discovered a background which may prevent this process from ever fulfilling 
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its potential in the study of CP-violation, i.e. Ki e~^e 77. 

AGS-845 also made major contributions to the study of — > e'^e~7 and Kl 
e^e'e^e" . This experiment is now completed. Its results are summarized in Table ^ 



Table 5: Results of AGS-845 



Mode 


Result 


Comments 


Kl 7r°e+e- 


< 5.5 X 10-9 


search for new scalars 
non-S.M. CP-violation 
future: S.M. CP-violation 


Kl e~^e~j 


(9.1±0.4+j;;J) X 10-6 
= -0.28 ± 0.083 


c.f. (9.6 ± 0.4) X 10-6 (theory) 
i.e. something besides p needed 


Kl e+e^77 


(6.6 ±3.2) X 10-' 116 


c.f. 5.8 X 10-^ for k* > 5MeV 
background to Kl — > vr°ee 


Kl e~^e~e'^e~ 


(3.04 ± 1.24 ±0.26) x 10"^ 85| 


c.f. 3.6 X 10-s 
background to Kl e~^e~ 



7.2 KEK-162 



KEK-162, shown in Fig. has been built by a KEK-Kyoto collaboration to pursue 
Kl — > iT^e~^e~ to ~ 10^^*^ ||118|| . It is quite similar in concept to AGS-845, with 
adaptations to the lower beam energy such as a more compressed layout, nitrogen 
rather than hydrogen in the Cerenkovs, etc. The main innovation is the electromag- 
netic calorimeter constructed of undoped Csl, a fast, bright scintillating crystal with 
excellent resolution. This calorimeter is designed to achieve 2% rms at 1 GeV. Most 
of the light in pure Csl is emitted with decay time < 30 nsec. Fast time response 
is crucial since the proponents of KEK-162 plan to expose their detector to an order 
of magnitude higher Kl flux than was seen by AGS-845. The daunting rates also 
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motivate the design of drift chambers which feature small cells, fast gas and custom 
TDCs. 

This experiment is currently setting up. 
7.3 FN AL- 731/799 

FNAL-731 was originally built to measure e'/e in K'^ 27r decays. In spite of being 
highly optimized for this purpose, it has produced several estimable rare decay re- 
sults. The apparatus, built by a Chicago/Elmhurst/FNAL/Princeton/Saclay group, 
is shown in Fig. Two nearly parallel Kl beams entered a 37m long evacuated 
decay region. A regenerator was shuttled between the beams on a pulse to 
pulse basis to provide Ks decays. A plane of thin trigger scintillators was situated 
approximately halfway down the vacuum decay vessel. Following the downstream 
decay region was a 4-station drift chamber dipole spectrometer, additional trigger 
hodoscope planes and an 804-element lead glass Cerenkov array. Downstream of this 
were photon vetoes, a 3m thick steel muon filter and, finally, muon veto hodoscopes. 
An extensive photon veto system bordered the acceptance. 

A comparison with detectors designed for lower energy beams is instructive. The 
FNAL-731 spectrometer was relatively more compact than the typical BNL or KEK 
detector and so had larger acceptance at the cost of worse charged track momentum 
resolution. Conversely, the resolution for photons was better for the higher energy 
experiment because of the 1/\/E behavior of the stochastic resolution term in lead 
glass. In FNAL-731, the decay region was shorter relative to the mean Ki decay 
length, diluting to some extent the advantage in acceptance. One major advantage 
of FNAL-731 was the clean beam in which the neutrons were of the same order as 
the Kl instead of ten or fifty times more numerous. This kept the detector rates 
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relatively low. 

A number of rare decay results were obtained from the seven-month 1987-88 run. 
Subsequently, the detector was reconfigured for the first stage of FNAL-799, a dedi- 
cated rare decay experiment. These changes included the removal of the regenerator, 
an upstream absorber, and the trigger scintillator plane in the vacuum decay vessel. 
These resulted in a large gain in sensitivity/ incident proton. The number of inci- 
dent protons/spill was also increased, so that the overall sensitivity/spill increased 
by more than an order of magnitude. Other modifications were an additional muon 
hodoscope plane and the development of an online processor which allowed the use 
of drift chamber information in the second level trigger. For a small portion of this 
run, a pre-shower detector was installed upstream of the lead glass array in order 
to improve the sensitivity to Kl 7r°77. The collaboration was also modified to 
consist of Chicago, Elmhurst, FNAL, Illinois, Colorado, UCLA, Rutgers, and Osaka. 
The experiment ran for 10 weeks in late 1991 and early 1992. 

The results of FNAL-73 1/799 on rare K decays thus far are summarized in Table 
The last three (preliminary) results come from the first stage of FNAL-799. This run 
is also expected to yield new results on — * 7r°e"*"e~, n^uV, and vr°77, 

as well as the first results of this program on Kl e"^e~7, Kl — > TTTree, Kl /z/xee, 
Kl — > ee77, 7r°/ie and other rare decays. 

The second stage of FNAL-799 is expected to begin taking data in 1994 or 1995. 
The experiment will be moved and a new beam line will be built which should allow 
higher primary intensity. The lead glass array will be replaced with one consisting 
of undoped Csl. The trigger planes and photon vetoes will also be upgraded. With 
these improvements, the experiment expects to obtain sensitivities < 10~^° for many 
rare decays. 
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Table 6: Results of FNAL-731/799 



TV T„ ^ „ 

iVlocle 


Result 


Comments 


Kl ^ IT e^e 


< 7.5 X 10-9 [ 


I14| 


searcn lor new scaiars, or viol. 


„o ^+ ^— 


< 4.5 X 10-5 [ 


119] 


needed for interpreting Kl — > vr'^ee 


Kl — ^ TT vr 7 


(3.19 ±0.16) X 10-5 [0 


direct emission, > 20MeV 


i^i, — > 71*^77 


(1.86 ±0.60 ±0.60) X 10-6 |7l| 


spectrum agrees with xPT 




< 2.2 X 10-^ 


HI 


CP-violating 


Kl ^ 7rV+/^- 


< 1.7 X 10-^ Q 


preliminary 


i^"/, e+e^e+e^ 


(4.17 ±0.83) X 10-* [0 


28 events 




(3.88 ±0.32) X 10-' [g 


167 events 



7.4 NA31 



NA31, although taking an approach to measuring e'/e which is radically different from 
that of FNAL-731, has had similar success in the pursuit of rare K decay modes. The 
detector, shown in Fig. |2^ was built by a collaboration of CERN, Dortmund, Edin- 
burgh, Mainz, Orsay, Pisa, and Seigen. The most striking difference of this detector 
from the others described in this section is the absence of a magnet. Charged par- 
ticle trajectories are determined by two planes of drift chambers, but their energies 
are measured by electromagnetic and hadronic calorimeters. The electromagnetic 
calorimeter, which is of the lead/liquid Argon type, has extremely good energy and 
position resolution. It is finely granulated and can distinguish two photons from one if 
they are separated by more than 1 cm. This configuration results in good acceptance 
for many rare decays, good energy resolution for electromagnetic particles, but rela- 
tively poor energy resolution for pions (e.g. the fractional resolution on a 50 GeV/c 
pion is ~ 9%). Other disadvantages of this configuration, such as the difficulty of 
distinguishing final state particles which are close in direction, are somewhat miti- 
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gated by the fine segmentation of the electromagnetic calorimeter. Similarly the lack 
of an energy/momentum comparison for distinguishing electrons from heavier parti- 
cles was partially alleviated by the longitudinal segmentation of the electromagnetic 
calorimeter and the by presence of the hadronic calorimeter. In addition there were 
dedicated particle identification systems such as a transient radiation detector and a 
muon filter. Completing the apparatus were triggering hodoscopes and photon veto 
counters. 

NA31 took data in 1986, 1988, and 1989. The rare decay results it obtained are 
summarized in Table 0. 



Ta 


ble 7: Results of CERN NA31 


Mode 


Result 


Comments 




< 4 X 10-^ [113] 


search for new scalars 


Kl —>■ 71*^77 


(1.7 ±0.3) X 10~^ 1 


70 1 


spectrum agrees with xPT 


Ks 77 


(2.4 ± 1.2) X 10-6 


38 


rate agrees with xPT 


Kl — >■ e"'"e~7 


(9.2 ±0.5 ±0.5) X 10"*^ |79| 
aK = -0.28 ±0.13 


cf. theory at (9.1 - 9.5) x 10^^ 
i.e. something beyond p needed 


Kl —>■ e+e^e+e^ 


(4 ±3) X 10'^ |84 


1 


First observation, 2 events 


K+ TT+XO; XO ^ e+e- 


< 6 X 10-^ to < 10-^ 


|120| 


limits depend on mx and Tx 



NA48, a successor to NA31, is now under construction. It features a fast liquid 
Krypton-based electromagnetic calorimeter and a magnetic spectrometer for charged 
particle momentum measurement. It should have an order of magnitude better sen- 
sitivity than NA31 to most rare decays. 
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7.5 AGS-777/851 

Fig. p4| shows the apparatus employed by a BNL/PSI/Washington/Yale group in AGS 
Experiments 777 and 851. The primary objects of these were respectively a search 
for 7r+/i^e~ and a study of —>■ 7r+e+e~. Other processes sought or studied 

in these experiments were tc^ /^^e^, 7r° —>■ e+e~, and e^e^ (where is a 

new light particle). 

A 6 GeV/c positive beam containing about 10^ / AGS pulse impinged on a 5m 
evacuated tank, wherein about 10% of the K'^ decayed. As the beam was unseparated, 
constituted only ~ 5% of the flux, the majority consisting of a roughly equal 
mixture of tt~^ and protons. The first element of the detector was a dipole run at a pt 
kick of 155 MeV/c. This served to remove the daughter products from the hot beam 
region and to separate them according to their charge. The beam then passed through 
holes and deadened regions in subsequent detector elements. Downstream of the first 
dipole was an MWPC spectrometer with four measuring stations. Two preceded and 
two followed a dipole run at Apx ~ 150 MeV/c with sense opposite that of the 
first dipole. Situated between each pair of measuring planes was an atmospheric gas 
Cerenkov counter. Downstream of the spectrometer were triggering hodoscopes, a 
lead-scintillator sandwich electromagnetic shower detector, and an iron/proportional 
tube chamber muon identifier. 

Since positive muons only were sought, the muon identifier needed to cover only 
the right (+) side. This was one of several optimizations made possible by confining 
the experiment to the n'^fi'^e' charge combination. On the left, where electron purity 
was more important than efficiency, H2 was used as C counter gas; whereas on on the 
right, where positrons had to be efficiently vetoed, CO2 was used instead. The choice 
of final state was also of great benefit to the trigger, since e~ are far less common in 
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decay than are (the most copious source of in decay is TT^e^u with 

branching ratio 0.0482; the most copious source of e~ is —>■ 7r+7r°;7r° — * e+e~7 
with product branching ratio 2.54 x 10"'^). 

The most dangerous backgrounds to 7r+/i+e^ stem from tt+tt+tt^ 

and Tr'^vr^; 7r° 'ye~^e~ decays. In the former, this can occur through various 

combinations of pion decay and misidentification. There are also several ways in 
which the latter process can mimic vr+yU+e^, the worst being the case where 

the TT"*" is mistaken for a /i"*", and the for a vr^. Since in this instance, a pion mass is 
misattributed to the e"*", the measured 3-body effective mass can exceed Mk^. Thus 
the usefulness of kinematic rejection is limited, so that powerful particle identification 
techniques are required. 

— > TT+TT+TT^ and tt+tt'^; tt^ 7e+e^ decays were not totally inimical to 

this experiment: they also served to calibrate and normalize it. The former process 
was used to evaluate the performance of the particle identification systems, design 
the geometrical reconstruction and kinematic fitting procedures, etc. 

For AGS-851, the C gas on the right side was changed from CO2 to The 
7r+yU+e~ trigger was dropped, and the requirement that the 7r+ be detected on 
the right was lifted. 

This program finished taking data in 1989. The AGS- 777 data analysis is now 
complete, that of AGS-851 continues. The results of these experiments thus far on 
rare K decay are summarized in Table |[ 

Subsequent to the completion of this program, the AGS Booster came on line, and 
the prospect of much greater available flux motivated the proposal of a successor 
experiment, AGS-865. The institutions collaborating are BNL, INR-Moscow, Dubna, 
New Mexico, PSI, Basel, Pittsburgh, Tbilisi, Yale, and Zurich. Since AGS-777 was 
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Table 8: Results of AGS- 777/851 



Mode 


Result 


Comments 


r^-f- -1- — -t- 
— s> TT^e jJ,^ 


< 2.1 X 10 ' |52 


Mh > 57 TeV 




(2.75 ± 0.23 ± 0.13) x 10"^ [|105| 
A = 0.105 ±0.035 ±0.015 


500 events 
suggests Ks — > n^e~^e~ small 


K+ ^ TT+XO; XO ^ e+e- 


< 1.1 X 10"^ |105| 

< 4.5 X 10-^ |121 


150 < rriee < 340MeV 
lOOMeV < ruee 
for Tx < 10^^^ sec. 



limited primarily by beam-associated background, a new beam has been designed to 
yield seven times more with no greater random rates than those of its predecessor. 
The detector, shown in Fig. |2^, is very similar to that of AGS-777. The geometrical 
acceptance has been increased by about a factor three with respect to that of the 
earlier detector, however, and the muon identifier covers both sides of the apparatus. 
Improvements to the background rejection power of the experiment include a fourth 
PWC plane at each measuring station, the use of aluminum HV wires to reduce mul- 
tiple scattering, an upgraded electromagnetic calorimeter, finer longitudinal sampling 
in the muon identifier, etc. 

The increases in the flux, the geometric acceptance, the triggering and recon- 
struction efficiencies and in running time are expected to yield a factor 70 improve- 
ment in the sensitivity of AGS-865 over that of AGS-777/851. This would allow a 
—>■ n^fi'^e sensitivity of 1.3 x 10~^^/event for example. At the same time samples 
of tens of thousands of decays such as vr+e+e^, 7r+/i+/i^ and 7r+77 should 

be accumulated. There are also a number of other interesting processes which could 
be studied with special runs and/or modest upgrades to the detector. These include 
CP- violating asymmetries in — > tt'^tt+tt^, T- violating polarization in K% decay. 
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and parity-violating polarization asymmetry in K'^ — > ir'^fj,'^/! . 
7.6 KEK-137 

KEK-137, a Tohoku/Tokyo/Kyoto/KEK collaboration, is one of the two most recent 
Kl — >■ 2-lepton experiments. These experiments have had to meet very significant 
challenges to achieve sensitivities significantly better than 10~^''/event. The fraction 
of Kl that decay and can be accepted in a practical sized apparatus is typically 
~ 1 — 2 X 10~^. To get enough Kl, the experiments have to work in the forward 
or near-forward direction where the neutron fiux is 10 — 100 times higher than the 
Kl. This translates into neutral beam fluxes of order 10^ per spill. Chamber and 
trigger plane rates are typically many MHz, and yet they must perform extremely well 
because to reject background, searches for Kl — > /xe must have excellent kinematic 
resolution. Since the primary background, K^^ decay followed by tt — > // via decay or 
misidentiflcation, occurs at the few % level and is topologically identical to the signal, 
particle identification power is of limited value. Most of the background rejection 
comes from kinematic and geometrical cuts. One exploits the fact that in the absence 
of measuring errors, the e "/x" pairs have effective masses less than Mk — 8.4 MeV. 
Both recent experiments kinematically and geometrically over-constrain their events. 
The momentum of each track is measured twice spectrometrically and in the case of 
muon candidates, once more via range. 

Here again, the background decays {K^s and for the Kl — > A*"^/^" case, K^^) are 
not all bad. In this case they serve to calibrate the particle identification devices. 
The relatively copious Kl tt+tt" are used to normalize the experiments as well as 
to calibrate the spectrometers. 

The double-arm spectrometer used by KEK-137 to search for Kl ^e, Kl 
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e^e~ , Kl l^^l^^, and other rare decay modes is shown in Fig. ^ A beam of 
~ 10'' Ki/pulse was made by directing a 1 — 2 x 10^^ proton beam from the KEK 
PS onto a 12 cm-long Cu target. The beam passed through a number of colhmators 
and sweeping magnets into a lOm-long evacuated decay volume. About 8% of the 
Kl between 2 and 8 GeV decayed in this volume. The rest of the neutral beam, 
which included about 10^ neutrons and 7s per spill, was conducted in vacuum be- 
tween the two spectrometer arms. The pt kick on each arm was 238 MeV/c, divided 
equally between the two dipoles. Daughter tracks in the Jacobean peak of the desired 
two-body reactions were consequently bent approximately parallel to the arm axes. 
Imposing a parallelism requirement greatly reduced the relative number of three-body 
decays accepted by the trigger. The resolution of this spectrometer for the calibra- 
tion Kl —>■ TT+TT" decays was 1.3 MeV/c^. As mentioned above, there were separate 
momentum measurements by the front and rear sections of the spectrometer and a 
third, coarse, measurement via a muon range array. The multiple measurements were 
primarily aimed at rejecting pions which decayed to muons in the detector. These 
were more dangerous than punch-through pions since the decay could disturb the mo- 
mentum measurement as well as the particle identification. Electrons were identified 
via atmospheric Cerenkovs filled with atmosphere, and with planes of lead-scintillator 
shower counters. 

This experiment ran for a period of about 2^ years, finishing data-taking in May, 
1990. The results are summarized in Table 0. 



7.7 AGS-791 

The second of the two-lepton experiments was AGS-791, a collaboration of UC- Irvine, 
UCLA, LANL, U of Pennsylvania, Stanford, Temple, and William & Mary. The de- 
tector, shown in Fig. had many similarities to that of KEK-137: a double arm 
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Table 9: Results of KEK-137 



Mode 


Result 


Comments 


Kl fie 


9.7 X 10"^^ 41| 




Kl e+e- 


1.6 X 10-1° |41 


one event in signal region 


Kl /i'^/u" 


(7.9 ±0.6 ±0.3) X 10-9 [|100U 


178 events 


Kl e+e~e+e" 


(6±2± 1) X 10-^ m 


partial reconstruction 



- double measuring spectrometer, electron ID via atmospheric C counter plus elec- 
tromagnetic shower counter, muon range array, etc. However there were significant 
difi^erences. The primary beam was about twice that of the KEK accelerator, and the 
daughter tracks entering the apparatus were roughly twice as stiff as those of KEK- 
137. Thus the muon identifier was considerably thicker and the Cerenkov needed to 
be filled with a gas of lower index of refraction (He-Ne mixture vs air). In E791, the 
arms shared their spectrometer magnets, and these were set to opposite polarities 
(the pt kick of each was ~ 300 MeV/c). Other differences from KEK-137 included 
a shorter decay volume but larger geometrical acceptance, lead glass instead of lead- 
scintillator sandwich counters, finer sampling in the muon identifier, etc. The Kl 
fiux impinging on the detector (~ 5 x 10^/ 1-second spill) was the highest yet used 
in a Kl decay experiment, and put severe requirements on the triggering and data 
acquisition systems. 

The analysis fiow was also similar to that of their KEK competitors. Normaliza- 
tion and spectrometer calibration was done with Kl tt+tt- decays, calibration of 
particle identification devices via Ki^ decays. 

In the course of analyzing the experiment, an unanticipated potential background 
to Kl — > fie was discovered. This was Kl vr^e^z/ decay in which both charged 
daughters were misidentified (vr as e, e as /i). Since a much higher mass is attributed 
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to the electron, the reconstructed two-body effective mass can span the region of the 
signaL Fortunately the particle identification power of the experiment was sufficient 
to reduce this background to < 10"^^. 

The data was taken in three runs over the period 1988-1990. The sensitivities 
reached were the highest ever attained in a decay experiment. The results are 



summarized in Table 10 



Table 10: Resu 


Its of E791 


Mode 


Result 


Comments 


Kl fie 


3.3 X 10-" 42 




Most sensitive K experiment yet 


Kl e+e- 


4.7 X 10-" 42 






Kl fi^fi" 


(7.0 ±0.5) X 10-9 


122 


718 events 



After the completion of data taking on AGS-791, a collaboration of UC- Irvine, 
Stanford, Temple, Texas, and William & Mary proposed a successor experiment, 
AGS-871. A schematic of their detector is show in Fig. |2^. Although many of the 
elements of the previous experiment will be reused, there are important differences 
in the design. The most striking of these is that instead of allowing the beam to 
pass unimpeded between the arms of the detector as in AGS-791, here it is stopped 
by a plug in the first spectrometer magnet. This allows larger geometric acceptance 
and lower rates in the downstream chambers and particle identification devices, at 
the expense of some increase in rates in the chambers near the plug. To help cope 
with this, the forward drift chambers will be replaced by high-rate straw trackers. 
A second important change is an adjustment in the spectrometer magnet fields so 
that there is a net Ap-r of 220MeV/c, to bend the two-body decay daughter tracks 
parallel to the beam, allowing a faster, more effective first level trigger and simpler 
Cerenkov optics. The over-bend (Ap^- = 440 - 220MeV/c) also leads to about a 
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20% improvement in the two-body effective mass resolution. Other changes include a 
longer decay volume, increased chamber redundancy, the use of H2 in the Cerenkov, 
better muon range resolution, and upgraded triggering and data acquisition systems. 
The Kl beam line will be lengthened to allow improved coUimation. 

Assuming a four-fold increase in the available AGS intensity, the overall improve- 
ment in sensitivity expected is about a factor 20. This implies single event sensitivities 
better than 10~^^. Thus a sample of over 10, 000 Kl will be accumulated as 

well few examples of — e^e . 

7.8 AGS-787 

The apparatus built by a BNL/Princeton/TRIUMF collaboration to carry out the 
first stage of AGS-787 is shown in Fig. |2^. The solenoidal configuration, unique among 
the detectors described in this section, was mandated by the problematic signature 
of n^uV. A 7r+ is hardly a novelty in the final state of a decay and 

since it alone of the three daughters is detectable, one has only the weak kinematic 
constraint < {Mj. — M^)/2Mk- What makes the experiment possible is the good 
signature of the backgrounds. The leading decays by far are i^^2 and K^2, each of 
which features a single charged track of unique cm momentum (205 MeV/c and 236 
MeV/c respectively). Thus with good kinematic resolution one can reject the two- 
body backgrounds by a large factor. In addition one can veto on the photons from 
Kt^2 and on the identity of the muon in i^^2 • Other potential sources of background 
(e.g. K'^ — > i^^3, i^^3, etc.) are less copious and all have combinations of 

at least two of the three 'handles' mentioned above: kinematic signature, detectable 
extra tracks, and charged daughter 7^ 7r+. 

Like all previous searches for ir^uV, AGS-787 employs a stopping beam. 
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This allows direct access to the kinematic features of signal and background. Other 
advantages are the feasibility of large geometric acceptance and veto hermiticity, the 
powerful particle identification techniques available at low energy, and the very good 
ratio of useful decays to unwanted beam particles. The latter results from the 
pure separated beams which are practical at low energies and the fact that one can 
stop a relatively large fraction of a low energy beam. This is quite important in 
an experiment in which the signature is a single unaccompanied vr"*". 

In the design of AGS-787, great efforts were made to minimize the presence of 
"dead" material. Energy deposited or interactions undergone in such material can 
compromise the veto or confound the particle identification. For example, a tt^ whose 
scatter in the stopping target goes undetected can defeat the kinematic rejection of 
i^7r2, if the TT^ decay photons are also missed. 

The background due to such "down-shifted" Kt^2 decays led the experimenters 
initially to concentrate on the kinematic region with Pt^+ > 205 MeV/ c where the only 
sources of 7r+ more copious than the signal are are vr+e+e^ and Tc'^'-fy. 

These have branching ratios of ~ 3 x 10^^ and < 10"^ respectively. The rejection of 
electromagnetic particles is extremely good in AGS-787 (e.g. only 1 or 2 of 10^ vr" 
are missed), so that these do not constitute a significant problem. In fact in this 
region, the most difficult backgrounds have proved to be fi'^u with muons 

which interact or enter the dead regions of the detector before making their range. 

Thus far, AGS-787 has operated at stopping rates up to ~ 300,000/ beam 
spill. To obtain this rate, about 1.5M 800 MeV/c from the LESBl were directed 
onto a BeO degrader 54cm in length. The were accompanied by an approximately 
equal number of protons and by about 3M 7r+. The emerged from the degrader 
with about 100 MeV of energy and came to rest in a highly segmented scintillating 
fiber target ||123||. were required to decay at least 2 nsec after stopping to be 
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accepted. Charged daughters which leave the target within about 30° of the plane 
transverse to the beam were tracked a cylindrical drift chamber onto which a field 
of IT was imposed. These then entered a cylindrical array of plastic scintillation 
counters and PWCs (range stack) which served to measure their residual energy and 
range. The range stack was read out on both upstream and downstream ends so 
that the detected particles could be localized in three dimensions. The counters 
were instrumented with 500 MHz transient recorders which recorded all scintillation 
light over an interval of about 10 /isec. This allowed the characteristic tt fi ^ e 
sequence to be observed in the stopping counter. This is a powerful signature for vr"*", 
which was defeated in < 1 out of 10^ cases. Comparisons among range, energy, and 
momentum of the daughter track also constitute a very effective particle identification 
technique. Surrounding the range stack was a cylindrical array of lead-scintillator 
sandwich counters (barrel veto) which veto gammas in the central region. The gamma 
veto was completed by lead-scintillator sandwich end-cap counters. 

Calibration and normalization of the experiment was done primarily via and 
K^2. 



Table 11: Resuhs of E787 



Mode 


Result 


Comments 




< 1.7 X 10-9 |59| 


familon, etc. search; 


— s> 7r+z/z7 


< 5.2 X 10-9 |9l H 


constrains new physics 


— > Tr+zi"*"/!" 


< 2.3 X 10-^ |TUB| 


should be discovered in '89-'91 data 


— > /i+z//i+/x-. 


< 4.1 X 10-^ |10| 


Higgs hunting ground 


7r+77 


< 10-6 |75| 


should be discovered in '89-'91 data 


K+ TT+^O; X° 77 


< 10-^ m 


mxo < 90 MeV/c2 



The results of AGS-787 are summarized in Table 



11 



These results are based on 
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data collected in an 1988 engineering run and in the first major physics run (1989). 
There were runs of about equal sensitivity in 1990 and 1991. During the 1989 run it 
was found that the instrumentation of the experiment was sufficient to reject back- 
ground in the kinematic region below the ir^2 as well as above it. Thus far the 
acceptance in the two kinematic regions has proved to be about equal. 

In 1989 a major upgrade of the beam and detector was approved with the aim of 
reaching the 10^^*^ sensitivity level necessary to probe the Standard Model predictions 
for TT^fV. The collaboration was augmented by groups from INS/Tokyo and 

KEK. A new low energy separated beam (LESB3) which provides much improved 
flux and purity was constructed. This beam was designed to exploit the AGS upgrade 
to deliver up to 1.5 x 10^ 800MeV/c X+ with a X+/7r+ ratio of 2/1. Extensive 
detector improvements were also undertaken. These include upgrades to the 
and TT"*" measuring devices, the photon vetoes, the electronics and data acquisition 
system. A new, brighter, stopping target was constructed out of close-packed 5mm 
square cross-section scintillating fibers. A new central drift chamber with only | 
the mass of its predecessor is being constructed. In the range stack, the scintillator 
read out granularity is being increased and the embedded proportional chambers are 
being replaced by far less massive straw chambers. Pure Csl end cap vetoes are being 
built to replace the current lead-scintillator sandwich devices. A pure Csl finer will 
be installed inside the current barrel veto. A number of supplementary vetoes will 
increase the hermiticity of the the detector. Transient recorders will be installed on 
the target and on all vetoes. Finally, a new trigger and data acquisition system, 
capable of taking ten times the previous rate is being developed. 
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8 Conclusions 

The study of rare kaon decays continues to be a crucial arena for the testing of 
electroweak theory. The decays K — > TriyV, can be computed rehably, and their study 
will yield valuable information on the CKM parameters p and rj, that will eventually 
permit us to test the three generation structure of the standard model. Other rare 
decays that are sensitive to these parameters are Kl 7r^£~^£~., —>■ vr+f+f" 
and perhaps Kl — > /^"''/^~ as well. Of course, these modes also constitute ideal 
candidates to search for new physics in the form of deviations from the standard 
model expectations. This is particularly true for CP violation. 

The search for the forbidden lepton flavor violating decays constitutes one of the 
simplest and most cost-effective ways to constrain interactions beyond the minimal 
standard model. The present level of sensitivity of these experiments is already testing 
energy scales that we will not be able to probe directly for many years. 

In the process of searching for the very rare and forbidden decay modes of the 
Kl and the K^, large samples of other, less rare, decays are being collected. Among 
them are the radiative decays that are dominated by long distance contributions. 
Their detailed understanding is crucial in the effort to use modes hke Kl vr^e+e" 
or Kl IJ''^IJ'~ to measure the parameters p and rj. These radiative decay modes 
are also interesting in their own right, in that they allow us to test the framework 
of chiral perturbation theory, and in that context they provide information on the 
non-perturbative aspects of the strong interactions. 

There are a number of accelerator developments under way to augment the supply 
of kaons for these studies. The AGS upgrade is well along. It will provide a fourfold 
increase in what is already the world's most intense source of kaons. Construction 
has begun on the Fermilab Main Injector which promises to provide a source of 
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comparable intensity at higher energy. Construction of the DA$NE storage ring at 
Frascati is also under way. This facility will provide a somewhat less intense source 
of kaons, but one in which the initial state of the kaons can be tightly controlled. In 
the somewhat more distant future, facilities such as TRIUMF Laboratory's proposed 
KAON complex can provide further large increments in sensitivity. 

The field of rare kaon decays has provided many discoveries of the highest im- 
portance. From an historical perspective, the reach in sensitivity of the present and 
near-future experiments is quite large. It would be very surprising if further discov- 
eries did not await these initiatives. 
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Figure Captions 

Fig. 1: Results from recent searches for lepton flavor violation in K decay, a) 
m(yue) vs ^^fjneup for Kl lie candidates from KEK-137; b) m{fie) vs 
for Kl — > fie candidates from AGS-791; c) mij^jjie) vs vertex quality 
variable for vr^/i+e^ candidates from AGS-777. 

Fig. 2: Box diagrams giving rise to — > ix^e"^ and K vr/i^e^ in models with 
massive neutrinos. 

Fig. 3: Short distance contributions to K ttuu. The full circle represents the 
effective one-loop sdZ coupling. 

Fig. 4: Potential long distance contributions to n^uV. 

Fig. 5: Search regions for vr+z/P from AGS 787. a) > Pktt2 and b) 

P7T+ < PK7T2- 

Fig. 6: Loop diagrams that give rise to 7r°77 at order The same 

diagrams without the line give rise to Ks — > 77. The X represents a 
vertex from Eq. |2.9|. 



Fig. 7: Examples of 0{p^) contributions to Kl 7r°77 in VMD models, a) Pole 
diagrams and b) Direct weak counter-terms. In both cases the X represents 
a weak transition, but in (b) it is 0{p^). 

Fig. 8: Rate for K^ — > 7r°77. a) Theoretical spectra. The solid line is the 0{p'^) 
xPT (ay = 0) result, and the dashed and dotted lines show the range for the 
values of ay given in the text. Phase space and pure VMD spectra are also 
shown; b) Data from Ref. [|70[ (solid histogram), compared with O(p^) result 
(dotted histogram). The latter has been multiplied by the experimental 
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acceptance (shown as crosses). Dashed histogram is calculated background. 

Fig. 9: Pole diagrams that give rise to C{z, v) in —>■ 7r^77. 
Fig. 10: Pole diagrams contributing to 77 

Fig. 11: Pole diagrams contributing to the direct emission in — > 7r+7r^7. 

Fig. 12: Short distance diagrams giving rise to and K'^ 7T~^i'^i~ . The 

full circle represents the effective one-loop sdZ vertex. 

Fig. 13: Dominant long distance contribution to . The vertical dashed 

line represents the cut to obtain the absorptive part. 

Fig. 14: Diagrams contributing to K Tc-f* at The full circle represents a 

vertex from Eq. p.7| , whereas the full box represents a vertex from Eq. p. 12 
The X is a 0{p'^) weak transition from Eq. p.9| . 



Fig. 15: Calculated decay distribution for ir^e^e , for values of men- 

tioned in the text. 

Fig. 16: Fit of —>■ 7T~^e~^e~ spectrum and branching ratio to predictions of xPT. 



From Ref. 105 



Fig. 17: Short distance contributions to the direct CP violation in Kl — > n^e^e 
Again, the full circle represents effective one-loop couplings. 

Fig. 18: Two photon contribution to the CP conserving absorptive part of - 

Fig. 19: Background to the process Kl —* ii^e^e^ . 
Fig. 20: Plan view of AGS-845 detector 
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Fig. 21: Schematic of KEK-162 detector 

Fig. 22: Elevation view of FNAL-731 detector 

Fig. 23: Schematic layout of NA31 detector 

Fig. 24: Plan view of AGS-777 detector 

Fig. 25: Schematic of proposed AGS-865 detector 

Fig. 26: Plan view of KEK-137 detector 

Fig. 27: Plan view of AGS-791 detector 

Fig. 28: Schematic of proposed AGS-871 detector 

Fig. 29: Side elevation view of AGS-787 detector 



